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ABSTRACT 


This thesis is comprised of four independent computer 
programs and their related operating instructions. Each of 
these programs focuses on a particular facet of tactical 
missile design. The topics covered include guidance and 
trajectory calculations, rocket motor propulsion sizing, 
warhead design, and aerodynamic coefficient determination. 
The programs are developed from accepted mathematical proce- 
dures and are tailored to optimize operator interaction for 
educational purposes. This thesis is intended to be uti- 
lized as a supplement to the thesis Tactical Missile Concep- 
tual Desigg by D.R.Redmon, Naval Postgraduate School, Sep- 
tember 1980. 
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I. INTRODUCTION 


The programs contained in this thesis were assembled 
expressly to supplement the work of Dan Redmon in his the- 
sis, Jactical Missile Conceptual Design. Two of the 
programs, LPATH, the trajectory model, and LAERO1, the aéro- 
dynamic coefficients program, originated in Redmon's thesis 
and were expanded/modified for use on the Naval Postgraduate 
School's new IBM 370 computer system. The other two pro- 
grams were generated for this thesis and utilize the proce- 
dures and principles outlined by Redmon. 

The specific intention of these programs is to provide 
students of tactical missile engineering and design with a 
method of solving complex mathematical routines rapidly and 
interactively. Each of the programs request data which are 
likely to be used as design parameters for the topics con- 
cerned. The programs also allow repeated operation with 
input alteration capability, allowing the user to personally 
optimize his design. This approach was chosen to allow stu- 
dents to understand the relationships various input parame- 
ters have with the final solutions. 
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II. TRAJECTORY MODELS 


A. DESCRIPTION AND ORIGIN OF PROGRAM 

This program applies the principles of missile guidance 
laws to the terminal phase (the last 5 to 10 seconds) of a 
missile trajectory in order to determine the maximum normal 
acceleration on the missile for a given scenario. Of the 
three general guidance law categories, pursuit guidance is 
not included in the program capability. It has been found 
that pursuit guidance invariably produces a tail-chase situ- 
ation, greatly reducing an anti-air nmissile's effectiveness 
against maneuvering targets of Similar speed characteris- 
tics. Line-of-sight guidance and proportional navigation 
guidance are both options of the progran. 

Figure (II-1) shows a typical encounter geometry as 
reguired for this program. The encounter is considered to 
cccur entirely within a two- dimensional plane. No differ- 
entiation is required or assumed concerning the orientation 
of the encounter plane. The plane may be at any angle to 
the horizontal as desired by the program user. The refer- 
ence direction is an arbitrary choice by the program user. 
The angles shown are positive in value, however the progran 
does not require positive angles. If TAL were 150 degrees, 
it could also Le entered as -210 degrees. The IRA term 
Cepresents the initial range to the target. 

Tangential velocities (air speeds) of the missile and 
the target are considered by the program to be constant 
throughout the problem. Since the program concerns itself 
with only the final moments of a trajectory, this is a rea- 
sonable consideration. Target normal accelerations, when 
specified by the user, are also held constant throughout the 
prcblem for the same reason. The missile normal 
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ogram Variables 
line of sight angle) 

L (target alpha) 

L (missile alpha) 

TSP (target speed) 

MSP (missile speed) 

PRE (initaral range’ to: 
target from missile 
launch_point 

IRM (initial range to 

missile from 

launch point) 


Vm 


VA Angles are shown as 
positive quantities 


Reference direction 


Figure (II-1). Encounter geometry 
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acceleration, however, is a function of the guidance laws 
and will vary according to the resulting flight path 
reguirements. A constraint on the encounter is that the 
missile normal acceleration must be zero at the start of the 
prcoblena. 

The program "flies" both the missile and the target as 
siuple points in space with no consideration given to aero- 
dynamics. The missile will always strike the target dead 
center when it is given the proper speed advantage for the 
encounter since there is no provision for statistical miss 
analysis. The unbounded nature of the missile normal accel- 
eration allows the missile to turn as rapidly as necessary 
to hit the target. 

This program analyzes the given encounter by time incre- 
ment calculations. As is the case with all integrations 
conducted by incremental steps, the accuracy of the results 
is a function of the increment size. The results will 
become increasingly accurate as the time increment is made 
Smaller. However, as the time increment decreases in size, 
the length of the output becomes increasingly longer. The 
user must balance the desire for accuracy against the amount 
of time he wishes to spend on the computer terminal. 

The primary output is a tabular listing of the missile 
and target coordinates at each time increment. The coordi- 
nate frame is established within the encounter plane with 
the abscissa oriented along the reference direction. The 
problem stops once the missile has passed its closest point 
of approach to the target. Output then provides the time to 
iatercept from time of problem initiation and the maxinun 
acceleration the missile was required to endure. A Versatec 
plot of overlaying successive encounters is an optional 
output. 
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This program originated as two seperate BASIC programs 
written by Redmon (Ref. 1] for use on the HP 9830 desktop 
Calculator. It was subsequently translated into FORTRAN IV 
for use on the IBM 370 computer system. 


Be. USER INSTRUCTIONS 

If it is desired to abort the operation of this progran 
prematurely, two methods are available. When the program is 
waiting for data entry, press (ENTER]. When the program is 
not waiting for data entry but is processing, type "HX" and 
press [ENTER]. Both actions will return the terminal to CMS 
mode. 

When the screen becomes full, or "MORE...." appears in 
the status area, clear the screen by pressing {ALT] and 
Simultaneously. At several points in the prograag, 
the user will be directed to "CLEAR SCREEN AND INPUT O", 
This is to provide proper positioning of the output on the 
screen for ease of reading. If any other symbol than 
" cra" should appear in the lower left of the screen, 
press and continue. 

1. Turn the terminal on with the red Switch. 

2. When the large "NPS" appears after about 30 seconds, 
press [RESET], then press [ENTER]. 

3. When "CP READ" appears in the status area on the 
lower right of the screen, type "L nnnnP", where nonn is 
your 4-digit user number, then press [ENTER]. 

4. You will now be asked for your password. Type it in 
(the characters will not appear on the screen), then press 
[ENTER]. 

5. Your personal file must contain a PROFILE EXEC rout- 
ine with the appropriate Fortran GLOBAL statement. If it 
does not, type "GLOBAL TXTLIB FORTMOD2 MOD2EEH", then press 


[ENTER] . 
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6. Type wMLENK TOexxxxP 191 AS 192 RE"; where xxxx is 
the 4-digit user number for the project file, then press 
[Ewrea]. 

7. You will now be asked for the project file password. 
Type it in (the characters will not appear on the screen), 
then press . 


8. Type “ACCESS 192 B™" and press | ENTER]. 


9. Press {ALT] and simultaneously to clear 
screen. 

10. Type "LPATH" and press [ENTER]. 

11. Input the following data as it is requested on the 
screen by typing the data and then pressing [ENTER] . Ensure 
that the data is input as either decimal or integer as spe- 
cified. The terms in parenthesis below are the progran 
variable names. 


Parameter Units Value range 

Trajectory option (TITLE) none O=Line-of-sight 
1=Proportional 

navigation 

2=Both 

Time increment (DEL) seconds Larger than .0001 
times the problem 
time 

Navigation constant (NAV) none 225° to 655 

LOS Angle (LOS) degrees 0.0 to 360.0 

ritial target Eanes meters Larger than the 

from the launch site (IRT) misSile range 

Target speed (TSP) n/sec Larger than 0.0 

Target normal n/sec/sec Positive is to 

acceleration (TAC) target's ett 

pad tz al eesti ae to degrees 0.0 to 360.0 the 

line of sight (TAL) 

Missile speed (MSP) m/sec Larger than 0.0 

Initial missile Eoige meters Such that impact 

from launch site (IRM) occurs in less 


than 10 seconds 

12. If you desire to rerun the problem, or want to run a 
New problem, answer the questions appropriately when asked 
by the terminal after the execution of the current problem. 
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13. To receive the printout and plot of your encounters, 
answer "no" to rerunning or restarting the problem when 
asked by the terminal and follow the directions presented on 
the screen. 

14. Upon completion of the program, type "LOGOFF" and 


press : 
switch. 


15. Turn the terminal off with the red 


C. EXAMPLE PROBLEMS 


1. Example I7-A. lLine-of-sight Non-maneuvering Cross- 
ing Target 


Not to scale 


Launch CA 
point “ fo 


Figure (II-2). Non-maneuvering crossing target 


Rt=10000 meters 
Rm=9000 meters 
60=30.0 degrees 
Ato=-90.0 degrees 
at=0 

Vt=200 meters/second 
Vm=800 meters/second 
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OUTPUT OF EXAMPLE II-A 
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Table (II-1) is the corresponding computer output for 
this encounter. As indicated, the missile maximum normal 
acceleration is 

am= -32.00 m/sec/sec or -3.26 g's. 
Figure (II-3) is the Versatec plot of the engagement. 

2. Exapole IT-B. Pp ortion Navigation Maneuverin 

crossing Target 
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Figure (II-4). Maneuvering crossing target 


Rt=10000 meters 

km=9000 meters 

$0=30.0 degrees 

Kto=-90.0 degrees 

at=156.8 m/sec/sec (16 g's) 
Vt=200 meters/second 
Vm=800 meters/second 
Navigation constant=3.06 
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Table (II-2) is the computer output for this encounter. 
The missile maximum normal acceleration is 
am= 283.94 m/sec/sec or 28.95 g's. 
Figure (II-5) is the Versatec plot of the engagement. 
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D. PROCEDURAL FLOWCHART 


Do you want just. proportional 
navagation? 


Establish initial angular 
relationships and determine 
Ainitra dl sais 
X=Y coordinstes, for, missile 
and target 


Have you exceeded 9999 itera- 
tions? 


SEEeEaenenail een 


iar iraaiaas 


Print error message 


Calculate angular celation- 
ships between aS oe and--tas- 
je 


Calculate nornal acceleration 
on missile 


Determine range between mis- 
sile nod Carget 


Calculate new X-Y coordinates 
for missile and target based 
on speeds. ’ 
and angular relationships 


- 


oe - 


Je. lcas 
hee? ee 


2 


= 


Soh ke ae aoe 


Print output | 


No "po you want. just line-of- 
Sight? 
yes 


I END | 
Establish initial angular 
relationships and determine 

; nea X=. 
coordinates for missile Sear 


lb a ae 


Increment iteration counter 
Have you exceeded 9999 itera- 
{ions 

9 ee 
| Print error message 
i END | 
Calculate new X-Y positions 
for missile & target based on 
speed |. ; 
and angular relationships 
Determine range between mis- 
Sile & target 
| Increment time | 
Calculate missile normal 
acceleration requirment based 
On Change in . 
angular relafionships 
Calculate new angular rela- 
tionships 
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E. PROGRAM CHANGES 
1. Language Translation. 

The two BASIC pro rams contained in Reference 1 were 

translated into standard FORTRAN IV. 
2. Program Condensation. 

The two individua programs were combined to form a 
Single integrated routine which allows the user to choose 
either of the two guidance laws or both for a given encoun- 
ter. The two original projrams were meshed such that only 
the actual guidance law algorithms are separate routines, 
all input and output rout nes are now common. 

3) el Dput..and Awtput Facility. 

The data input instructions were modified to maxin- 
ize user facility on the I8M 370 System 3278 terminals. 
Completion of data input is now followed by data feedback 
for user verification prior to actual program execution. 

Data output has been expanded to provide data deliv- 
ery to three destinations; to the user's terminal for obser- 
vation, to a printer file for hardcopy duplication of the 
terminal display, and to a plot data file for subsequent use 
by the plot program. The iestinations are options chosen by 
the user for each execution of the program. Up to nine 
different problems can be printed and plotted each time the 
program is entered. 
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The program can now be rerun multiple times without 
exiting and re-entering each time. The user has a choice of 
either rerunning the same problem or initiating a new prob- 
lem completely. 

4. Plot Program. 

PATHPLOT FORTRAN was developed to produce a Versatec 
plot picture of the encounter. It will produce a single 
plot each time the program is entered and will plot up to 9 
engagements in an overlaying manner. This format was 
adopted to allcw comparisons of successive input data modi- 
fications. 

5. Data Qvercapacity Check. 

If the user initiates a problem requiring more than 
9999 iterations, the program will stop. The user will be 
notified of the error and given the opportunity to rerun the 
problem. 

6. Znitial Missile Acceleration. 

Initial missile acceleration was removed as a user 
input variable and established as zero. Due to the mathe- 
Matted! nature (i.e., nO physical constraints) of the pro- 
gram, any "wrong" accelerations of the missile in the ini- 
tial state were immediately "corrected" by the algorithn. 
The modei is better served by providing no initial accelera- 
=20ns. 

io ineta. Angles. 

Both the target and the missile theta angles (the 
angles between the velocity vectors and the reference direc- 
tion) were removed as user input variables. The program now 
Calculates the theta angles from other input data, reducing 
redundancy and possible contradiction of input data. 

38. [Initial Conditions Perspective. 

Originally, the missile guidance calculations 

Started at t=0. Specifically, the anchor point for the line 
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of sight solution was the point where the missile commenced 
the problem, whether or not that was on the launch site. 
This produced a situation removed from reality where the LOS 
anchor point should be at the fire control location, usually 
at or near the launch site. A similar, though less pro- 
nounced, condition existed for the proportional navigation 
sclution. The program was modified to provide proper posi- 
ticning of the external guidance reference. AS a result of 
the modification, an additional output is the correct lead 
angle for the missile at the start of the problem. This 
angle is based on the assumption that neither the target nor 
the missile have maneuvered prior to time t=0. 
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IIIT. WARHEAD DESIGN 


A. DESCRIPTION AND ORIGIN 

This program develops a warhead using the same methods 
as presented by D Redmon (Ref. 1]. However, its capability 
is somewhat greater and it applies the relationships in a 
Slightly different manner. This program was written for 
this thesis in FORTRAN IV for use on the IBM 370 computer 
system 3278 terminals. It is configured expressly for oper- 
ator-computer interaction. 

me Starts with an initial input of data from which a 
table of fragment initial velocities is generated and pre- 
sented to the user. From this, the user selects values for 
fragment mass and impact velocity and another table is gen- 
erated and displayed. From this one, the desired probabil- 
ity of hit given a detonation is selected and the final 
solution is produced. At various points during the opera- 
tion of the program, the user will have opportunities to 
alter or revise selected parameters. 

The program is limited to a cylindrical warhead with 
either a solid or hollow core. The fragments: are required 
to be square in shape and are sized by the program. Figure 
(IIIiI-1) illustrates the shape of the warhead and the loca- 
tion of various input and output quantities. 

Initially, the target altitude is used to determine the 
atmospheric density, temperature and speed of sound. These 
values are, in turn, used to calculate required initial vel- 
ccities for the fragments. The program is preloaded with 
Various values for the fragment mass and impact velocity, 
which are used to generate the table of initial velocities. 
These velocities come from the following relationship: 


V (hit) =Vif{exp (-ks) } 
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where V(hit) is the impact velocity, Vi is the initial vel- 
ocity, s is the kill radius, m is the fragment mass, Pa is 
the atmospheric density, A is the plan area of the fragment, 
and Cd is the drag coefficient of the fragment. 

The ballistic limit velocity is calculated for the vari- 
ous presized fragments and is provided as a reference when 
choosing an appropriate impact velocity. The ballistic 
limit velocity is that velocity at which one half of the 
fragments will penetrate the target's skin and the other 
half will not. The empirical relationship, developed by A. 
E. Fuhs (Ref. 6], presents the ballistic limit velocity as a 
function of the fragment size to skin thickness ratio. His 
function dealt with steel fragments impacting an aluminuao 
plate. His results were modified to qualitatively account 
for different skin materials and fragment densities. 

Next computed is the fragment spray angle and the criti- 
cal miss distance. The spray angle is a function of the 
initial velocity, the detonation velocity and the warhead 
length. The critical miss distance is defined as the range 
where the fragment spray exactly covers the entire target. 
The critical miss distance is used by the program to sepa- 
Tate the two functions which determine the average number of 
hits received by the target. The program assumes the target 
is always centered within the fragment spray and aligned 
perpendicular to the central ray of the spray. 

A selection of warheads is then produced, one for each 
of a preloaded set of Pd's (probability of a hit given a 
detonation) to provide the user with a Pd versus warhead 
weight/size trade off comparison. This sizing process is 
based entirely on the following relationship: 


Pd=1-exp (-a) 
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where a is the average number of hits. The variable a, as 
shcewn by Redmon (Ref. 1], is a function of the cube of the 
warhead radius. The user then chooses a desired Pd which, 
in turn, produces the final warhead sizing. 

Values that are shown in tabular form for user selection 
and input into the program are not limited to those dis- 
played. Any value in between the displayed values or con- 
pletely outside of the value range may be used. The one 
exception is that Pd can never be selected to be larger than 
-999 and may even, if forced by the program limitations, be 
reguired to be lower. Since ultimately in this algoritha, 
Pd is a function of the warhead radius, the maximum value 
for Pd may ke reduced in order to keep the radius within the 
Original input parameters. The user is notified if this 
condition occurs. 

Useful reference information for some common explosives 
and case materials is contained in the following tables. 


Table (III-1). Characteristics of common explosives 


ate aS set ae t/sec) Vd (ft/séc) 
RDX - 0596 9300. 26837. 
HMX - 0665 10230. 29934. 
FEIN - 0625 9300. 212316 
TETRYL «0585 8200. 25:7 55% 
COMP B - 0607 8800. 257122 
OCTOL - 0650 9500. 28356. 


Table (III-2). Densities of common case materials 


Alunjnun - 100 
Uranium -688 
Titanium 163 
Lead ~410 
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Figure (III-1). Porm of warhead as used by the program 
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Be. USER INSTRUCTIONS 

If it is desired to abort the operation of this progranm 
prematurely, two methods are available. When the program is 
waiting for data entry, press [ENTER]. When the program is 
mec ating for data entry but is processing, type “Hx" and 
press [ENTER]. Both actions will return the terminal to CaS 
mode. 

When the screen becomes full, or "MORE...." appears in 
the status area, clear the screen by pressing and 
simultaneously. At several points in the progran, 
the user will be directed to "CLEAR SCREEN AND INPUT 0". 
This is to provide proper positioning of the output on the 
screen for ease of reading. If any other symbol than 
"GEA" should appear in the lower left of the screen, 
press and continue. 

1. Turn the terminal on with the red Switch. 

2. When the large "NPS" appears after about 30 seconds, 
press [RESET], then press [ENTER]. 

3. When "CP READ" appears in the status area on the 
lower right of the screen, type "L nnnnoP", where nonnn is 
ycur 4-digit user number, then press [ENTER]. 

4. You will now be asked for your password. Type it in 
(the characters will not appear on the screen), then press 
(Eyres). 

5. Your personal file must contain a PROFILE EXEC rout- 
ine with the appropriate Fortran GLOBAL statement. If it 
does not, type "GLOBAL TXTLIB FORTMOD2 MOD2EEH", then press 
[enTER]. 

oe, Lype "LINK TO xxxxP 1917 AS 192 RR", where xxkxx is 
the 4-digit user number for the project file, then press 
(EwrEy] . 

7. You will now be asked for the project file password. 
Type it in (the characters will not appear on the screen), 


then press |ENTER]. 
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8. Type “ACCESS 192 BY" and press | ENTER]. 
9. Press and {CLEAR] simultaneously to clear 
screen. 


10. Type “LBOMB" and press [ENTER]. 

11. Input the following data as it is requested on the 
screen by typing the data and then pressing [ENTER] . Ensure 
that the data is input as either decimal or integer as spe- 
cified. The terms in parenthesis are the program variable 
nagwes. Input the following variables as decimal values: 


Explosive density (XDEN) iED/GUsL0 
Explosive Gurney constant (GC) ft/sec 
Explosive detonation velocity (VD) ft/sec 
Case material density (CDEN) Ep7euslt 
Desired kill radius (RKILL) feet 
Wachead diameter (DW) inches 
Warhead Llength-to-diameter ratio (LTD) 
Target length (LT) feet 
Target mean width (WT) feet 
Target vulnerability, P(k/h) (PKH) 
Target altitude (ALT) feet 
Target skin thickness (TTH) inches 
Target skin material (MAT) 1.0=Aluminun 
2.0=Fiberglass/Kevlar 
3.0=Steel 


12. After entering the above data, you will be presented 
with an initial velocity table built around your desired 
kill radius. The initial velocities will be listed as a 
function of fragment mass and impact velocity. Also pro- 
vided will be the ballistic limit velocities for each of the 
fragment masses. Input the following parameters as decimal 


values: 
Fragment mass (IFMLB(1)) grains 
Impact velocity (VHIT(1)) ft/sec 
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13. You will now be presented with a shopping list of 
warheads developed for a range of Pd's. The warheads are 
described by the following quantities: 


Warhead weight in pounds (total weight) 
Case thickness in inches 

Core diameter in inches 

Number of fragments from the warhead 
Number of fragments hitting the target 
Edge length of the fragments in inches 


Input the following parameter as a decimal value: 
Desired probability of a hit given a detonation (PDF) 


14. If you desire to rerun the problem, or want to run a 
new problem, answer the guestions appropriately when asked 
by the terminal after the execution of the current problem. 

15. To receive the printout of your encounters, answer 
"no" to rerunning or restarting the problem when asked by 
the terminal and follow the directions presented on the 
screen. 

16. Upon completion of the program, type "LOGOFF" and 


press ENTER]. 
i7. Turn the terminal off with the red Switch. 


C. EXAMPLE PROBLEMS 

tH. Exasple [III-A 

It was desired to build a warhead which would kill a 

typical cruise missile flying at 100 feet altitude. The 
warhead was required to have a lethal radius of 50 feet with 
a Pd of 0.98. TNT was selected for the explosive and steel 
was chosen for the case material. The diameter of the mis- 
Sile was 13.5 inches. 
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Table (III-3) outlines the input parameters. Takle 
(III-4) shows the program output. After the initial veloc- 
ity table was displayed at the terminal, the fragment mass 
was chosen to be 130 grains and the impact velocity was 
selected to be 2500 feet per second. After the Pd table was 
displayed, 0.98 was entered as the kill probability. 

2. ‘Example IiI=8 

A warhead was required which would kill a typical 
Manned aircraft at 30000 feet. A lethal radius of 75 feet 
was specified. The warhead was limited in weight to 50 
pounds and needed to have a core diameter of at least 5 
inches. The diameter of the missile was 10.0 inches. HMX 
was chosen as the explosive and depleted uranium as the case 
material. 

Table (III-5) outlines the input parameters. Table 
(III-6) shows the program output. After the initial veloc- 
ity table was displayed at the terminal, the fragment mass 
was chosen to ke 210 grains and the impact velocity was 
selected to be 5000 feet per second. After the Pd table was 
displayed, it was determined that the fragments were too 
large and the fragment mass was then reduced to 100 grains. 
The impact velocity was also reduced to 3000 feet per sec- 
cond. When the Pd table was redisplayed, 0.995 was chosen as 
the desired kill probability. 
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INPUT DATA FOR EXAMPLE III-A 
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OUTPUT DATA FOR EXAMPLE III-A 


TABLE (III-4). 
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INPUT DATA FOR EXAMPLE III-B 
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OUTPUT DATA FOR EXAMPLE III-B 


TABLE (III-6). 
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D. PROCEDURAL FLOWCHART 


Petites | 


earns: itt 


Calculate characteristics of 
atmosphere at altitude 


Calculate initial velocities 
for all combinations of preset 
impact 
velocities and fragment 
weights at the kill radius 


Print table of initial veloci- 
es 


- 


ae you want to change the kili ne 
radius 7 


YS 


Input desir i fragment mass 
and impact velocity 


ment spray radius 


Calculate the critical frag- | 


Es ehe kil radius larger than 
the Critical radiu 


Calculate average number of 
hats as a fanction of 1/7R 


Calculate average number of 
hits as a function of 1/8 
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si ae et Oe 


Calculate individual warhead | 
solutions for pease! values of 


trons 


If Pd was restricted due to 

max available warhead radius, 
Ene 

warning message 


Do you _want to change the 
fragment mass? 


eo @2e @ @e@e@ @ ee = oc3oe2 weer @ @w@e@ ewe @ @ ow 


Calculate refined value of 
charge-to-mass ratio 


Determine warhead solution for 
specified Pd 
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IV. PROPULSION MOTOR SIZING 


A. DESCRIPTION AND ORIGIN OF FROGRAM 

This program provides a method for the preliminary siz- 
ing of a solid propellant rocket motor for a boost-sustain 
trajectory of a tactical missile. The analytical method was 
developed by Redmon {Ref. 1] and was expanded with the addi- 
fron Of material from Platzek [Ref. 2] and Hill [Ref. 3]. 
The program was written for this thesis in FORTRAN IV for 
use on the IBM 370 computer system. Essentially, the analy- 
sis consists of sizing the booster motor from differential 
velocity and acceleration requirements with limitations 
imposed by the physical dimensions of the missile. The 
bocster is at all times considered to be a core-burning 
motor. The sustainer motor calculations are controlled ry 
the maximum range specified by the user and by the solution 
of the booster motor. The sustainer can be either a core- 
burning or an end-burning motor. 

The rocket motor configuration is assumed to be either 
an integral bocster-sustainer motor as shown in Figure 
(IV-1) or a staged motor as in Figure (IV-2). The booster 
and the sustainer always burn exclusively, or, in other 
words, one is not burning while the other one is. The 
nozzle half angle is specified by the user consistant with 
Space available in the missile. If a staged motor is uti- 
lized, both nozzles will have the same half angle. 

The booster calculations start by determining the amount 
of thrust needed to boost the total weight of the missile to 
its cruise velocity at the prescribed acceleration. From 
this, the necessary amount of propellant is estimated and 
the process is iterated to account for the decreasing mass 
Situation. The chamber pressure is initially estimated by 
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Figure (IV-1). Integral booster-sustainer motor 
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Figure (IV-2). 


Staged motor 
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minigizing the motor weight to propellant specific impulse 
ratio, as presented by Redmon (Ref. 1]. The next step of 
the program is to size the nozzle using the now known ini- 
tial pressure ratio. Once the ideal nozzle is developed, 
the chamber pressure is raised or lowered as necessary in an 
attempt to cause the nozzle exit diameter to exactly match 
the booster diameter. However, if the integral motor option 
is used, the program will drive the nozzle as small as it 
can without violating one of the following limits in order 
to increase the probability that the sustainer will operate 
properly. The iterative process hasS two limits: an abso- 
lute maximum of 2000 PSI chamber pressure [Ref. 2] and a 
Minimum of 1000 PSI if the pressure had previously been 
higher. The solution can be less than 1000 PSI if the pres- 
sure remained below that level throughout the problem. Also 
limited is the exit pressure to ambient pressure ratio. At 
the low end, it is limited to 0.4 to prevent flow separation 
wiethpe nozzle. At the high end, it is limited to 2.5 to 
prevent excessive underexpansion and loss of physical real- 
ity in the program resuits [Ref. 4]. The burn rate is 
initialized at a potential maximum of 1.25 inches/second 
(Ref. 2] and is allowed to decrease to arrive at a compati- 
ble burn area and web thickness combination. 

The sustainer motor, in the integral motor case, is vir- 
tually a continuation of the booster solution. The initial 
thrust requirement is determined by increasing the cruise 
Speed drag to account for speed loss through maneuvers. It 
is also then refined for weight change if any climbing or 
diving is required to reach the target altitude. The nozzle 
is the same one as developed for the booster except that the 
throat area will be expanded as a result of the erosion 
effect. From the new area ratio, a pressure ratio is deter- 
mined. The chamber pressure and thrust coefficient are then 
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iterated until a steady chamber pressure evolves to provide 
the required thrust. If at any time it drops below 125 PSI, 
the program will stop since this is considered a nininun 
chamber pressure for proper propellant combustion. The exit 
pressure to ambient pressure ratio remains subject to the 
same restrictions. The burn rate starts at 0.45 inches/sec- 
ond {Ref. 2] and is decreased to provide an acceptable web 
thickness and burn area. The solution can be either an end 
burner or a core burner, depending on the burn area 
required. 

The sustainer for the staged motor is solved in essen- 
tially the same manner as the booster. The two major excep- 
tions are that it can produce an end burner if the burn area 
is small enough and the thrust required is based on the 
cruise drag and not the velocity to be gained. The chamber 
pressure is limited to an absolute maximum of 800 PSI anda 
Minimum of 250 PSI if the pressure had previously been 
higher [Ref. 2]. The same pressure ratio restrictions apply 
and the burn rate starts at 0.45 inches/second. 

Other motor-nozzle combinations can be created from the 
results of this program on a first-order approximation 
basis. Figures (IV-3) and (IV-4) illustrate two possible 
methods for separating the motor nozzles without resorting 
to staging. Although these two methods will probably have 
the same weight disadvantage that staging does, they can be 
packaged in to a Smaller volume of Space. A hint for "crea- 
tive construction" is to rerun the problem after sizing the 
motors uSing a smaller missile diameter to force the nozzles 
to a smaller size. 
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Figure (IV-3). Separate nozzles (nonstaged motors) 
The values are obtained from the staged motor option. 
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Figure (IvV-4). Concentric nozzles (nonstaged motors) 
The values are obtained from the staged motor option. 
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B. USER INSTRUCTIONS 

If it is desired to abort the operation of this program 
prematurely, two methods are available. When the program is 
waiting for data entry, press fentTeR]. When the program is 
not waiting for data entry but is processing, type "HX" and 
press [ENTER]. Both actions will return the terminal to CMS 
mode. 

When the screen becomes full, or "MORE...." appears in 
the status area, clear the screen by pressing and 
Simultaneously. At several points in the progran, 
the user will be directed to "CLEAR SCREEN AND INPUT 0". 
This is to provide proper positioning of the output on the 
screen for ease of reading. If any other symbol than 
"GAR" should appear in the lower left of the screen, 
press and continue. A 

1. Turn the terminal-on with the red |O| switch. 

2. When the large "NPS" appears after about 30 seconds, 
press [RESET], then press [ENTER]. 

3. When "CP READ" appears in the status area on the 
lower right of the screen, type "L nannP", where onnn is 
your 4-digit user number, then press [enter]. 

4. You will now be asked for your password. Type it in 
(the characters will not appear on the screen), then press 

5. Your personal file must contain a PROFILE EXEC rout- 
ine with the appropriate Fortran GLOBAL statement. If it 
does not, type "GLOBAL TXTLIB FORTMOD2 MOD2EEH", then press 
[ENTER]. 

Seu ctype; "LINK. TO, xxxxP 191 AS 192 BR", where xxxx is 
the 4-digit user number for the project file, then press 
[Eyres] . 

7. You will now be asked for the project file password. 
Type it in (the characters will not appear on the screen), 


then press ° 
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8. Type "ACCESS 192 B" and press | ENTER|. 


9. Press j|ALT] and |CLEAR| simultaneously to clear 


screen. 


1. Lype “LPROP"® and press {ENTER . 


11. Input the following data as it is requested on the 


screen by typing the data and then pressing |ENTER]. 
that the data is input as either decimal or integer as spe- 


Ensure 
cified. The terms in parenthesis are the program variable 
names. 

The following variables are required inputs for both 
mctor option problems. 


Motor option (I MOTO R) He =integral motors 
1=staged motors 

Launch altitude (L ALT) feet 

Launch weight (WL) pounds 

Launch velocity (V BI) feet/second 

Launch elevation angle (ELB) degrees 

Boost acceleration (A) gravities 

Cruise velocity (V BF) feet/second 

Cruise velocity drag (DRAGS) pounds 

Maximum range to target (R) nautical miles 

Maximum target altitude (TALT) feet 

Bocster propellant 

specific impulse (I SP3) seconds 

Booster propellant ; 

density (DENSB) pounds/cu. inch 

Bocster exhaust. 

specific heat ratio (GB) 

Sustainer propellant 

Specific impulse (I SPS) seconds 

Sustainer propeilant ; 

density (DENSS) pounds/cu. inch 

Sustainer exhaust. 

specific heat ratio (G S) 

Nozzle half angle (ALN) degrees 
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The following variables are required inputs for the 


integral motors option only. 


Nozzle design altitude (ALTBN) feet 


Nozzle erosion rate (ER) inches/second 
Missile diameter (D B) inches 

Case yield strength (YIELD) PSE 

Case density (DENSC) pounds/cu. inch 


The following variables are required inputs for the 


staged motors option only. 


Booster design altitude (ALTBN) feet 


Booster diameter ( DB) inches 

Booster case 

yield strength (YIELD) PSI 

Booster case density (DENSC) pounds/cu. inch 
Sustainer design 

altitude ( ALTS N) feet 

Sustainer diameter (D) laches 
Sustainer case 

yield strength (YIELD) PSI 

Sustainer case density (DENSCS) pounds/cu. inch 


12. This program will cue the user when the input paran- 
eters have dictated a scenario which either cannot be 
achieved within reality or produce less than optimun 
requirements on the propulsion system of the missile. They 
are not definitive and are only intended to make the user 
aware of a situation which may need correction. The follow- 
ing is a list of available cue messages with short 
definitions. 

"SUSTAINER NOT CALCULATED SINCE THE BOOSTER BURNOUT 
RANGE EXCEEDS THE DESIGN RANGE." This can result fren, 
entering an extremely short range for the missile, or it 
Can be caused by a very low acceleration reguirement. 
“SUSTAINER NOT CALCULATED; THE BOOSTER NOZZLE DESIGN 
PREVENTS SUSTAINER OPERATION. RECOMMEND STAGING OR 
INDEPENDENT NOZZLES," This occurs only when eee the 
fo) 


integral motors option. Then scenario described the 
program can cause the booster nozzle to be too large to 


ey 


maintain the proper chamber pressures when the motor has 
shifted to sustainer operation. This usually occurs 
when a large acceleration is demanded but the thrust 
required fcr the cruise trajectory is small. 


“BOOSTER BURN RATE WAS ADJUSTED TO CORRESPOND WITH THE 
WEE THICKNESS.“ This occurs quite often and Simply 
indicates that the burn rate was decreased from 1tS 
petential physical maximum of 1.25 inches/second. 


"THE BOOSTER NOZZLE DESIGN WAS RESTRICTED DUE TO THE 
MISSILE DIAMETER," and "THE SUSTAINER NOZZLE DESIGN WAS 
RESTRICTED DUE TO THE MISSILE DIAMETER." fhe nozzle was 
not able to be designed for optimum pressure ratios at 
the mid point of the boost EeapeGe ony. Usually, the 
exit diameter is solved larger than the missile diameter 
and is subsequently reduced to fit. 


"BOOSTER NOZZLE DESIGN IS NOT OPTIMUM DUE TO EXCESSIVE 
BOOSTER CHAMBER PRESSURES," and "SUSTAINER NOZZLE DESIGN 
IS NOT OPTIMUM DUE TO EXCESSIVE SUSTAINER CHAMBER 
PRESSURES." If the nozzle cannot be downsized without 
aa cage e peceaure thresholds (2000 PSI for the booster 
and 800 PSI for the sustainer), the chamber pressure is 
held just relow the pressure threshold and the nozzle 
area ratio will be adjusted to allow the nozzle to fit 
in the missile. 


"THE SUSTAINER MOTOR HAS A CORE-BURNING GRAIN." The 

reguired burn area for the sustainer was too large tco 

peeart an end burning grain with a properly realistic 
urn rate. 


"SUSTAINER BURN RATE WAS ADJUSTED TO CORRESPOND WITH THE 
WEB THICKNESS." This indicates the sustainer burn rate 
was lowered from a potential maximum of 0.45 inches/sec- 
ond to provide a proper web thickness. 


"THE SUSTAINER MOTOR HAS AN END BURNING GRAIN." The 
required burn area for the sustainer was small enough to 
permit an end burning grain. fhe burn rate is then 
adjusted to correspond with the nonreduceable burn area. 


"REESTIMATION OF LAUNCH WEIGHT IS REQUIRED FOR THESE. 

MISSILE PERFORMANCE PARAMETERS." The scenarlo described 
to the program produced a motor whose weight is either 

eae han 75% of the total or less than 25% of the 
otal. 


"ENLARGEMENT OF DIAMETER IS RECOMMENDED DUE TO A VERY 
HIGH LENGTH-TO-DIAMETER RATIO FOR THE MOTOR." This cue 
indicates the length-to-diameter ratio is greater than 
15. Other components of the missile will make the mis- 
Ssile's overall length-to-diameter ratio even larger. 


13. Immediately after completion of the solution, the 


prcgram will ask if you want to receive a hardcopy printout 


of that particular solution. A "yes" answer stores that 


scluticn in a file for retrieval by the user when he fin- 


ishes with the progran. 
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14. If you desire to rerun the problem, or want to run a 
new problem, answer the questions appropriately when asked 
by the terminal after the execution of the current problem. 

15. To receive the printout of your encounters, answer 
"no" to rerunning or restarting the problem when asked by 
the terminal and follow the directions presented on the 
screen. 

16. Upon completion of the program, type "LOGOFF" and 


press [ENTER]. se 
Switch. 


i7.. -rurn the terminal off with the red 


(Ge EXAMPLE PROBLEMS 
1. Example IV-A. Integral motors, common nozzle 


The following parameters are input for the integral 
motor example; 


Launch altitude 35.0 feet 

Launch weight 1000.0 pounds 
Launch velocity 0.0 feet/second 
Launch angle 60.0 degrees 
Average acceleration 3:10 .0'cg*s 

Cruise velocity 4000.0 feet/second 


Drag at cruise velocity 1500.0 pounds 
Maximum range 20.0 miles 
Final (target) altitude 50000 feet 
Booster propellant 


specific impulse 260.0 seconds 
Bocster propellant ; 
density 02075 Lbs/cu. inch 
Booster exhaust specific 

heat ratio 1.244 
Sustainer propellant 

Specific impulse 210.0 seconds 
Sustainer propellant : 
density 0.065 Llbs/cu.inch 
Sustainer exhaust. 

specific heat ratio We27sU 

Nozzle half angle 20.0 degrees 
Nozzle design altitude 0.0 feet 
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OUTPUT OF EXAMPLE IV-A 


TABLE (IV-1). 
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Nozzle erosion rate 0.001 inches/second 


Missile diameter 10.0 inches 
Yield strength of 
case material 180000.0 PSI 
Density of case ; : 
material O.2662 iLbs/cu. inch 


The solution for this problem is presented in Table 


(IV-1). 
2. Example IV-B. Staged motors, separate nozzles 
The following are input for the staged motor prblen: 
Launch altitude 35.0 feet 
Launch weight 2200.0 pounds 
Launch velocity 0.0 feet/second 
Launch angle 30.0 degrees 
Average acceleration 25.0 g*s 
Cruise velocity 2200.0 feet/second 
Drag at cruise velocity 1000.0 pounds 
Maximum range 50.0 miles 
Final (target) altitude 75000 feet 
Booster propellant 
specific impulse 250.0 seconds 
Booster propellant ; 
density -0647 lbs/cu.inch 


Booster exhaust 


Specific heat ratio Ue.225 
ae ee ae 

specific impulse 205.0 seconds 
Sustainer propellant ; 
density ~0625-1bs/cu.inch 


Sustainer exhaust 


specific heat ratio 1258 

Nozzle half angle 15.0 degrees 
Booster design altitude 0.0 feet 
Booster diameter 14.5 inches 
Yield strength of 

booster case 180000.0 PSI 

Density of booster ; 
case material 022662 bs/cu..2anch 


Sustainer design altitude 
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OUTPUT OF EXAMPLE IV-B 


TABLE (1V-2). 
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Sustainer diameter 14.5 inches 


Yield strength of 
sustainer case 180000.0 PSI 


Density of sustainer ; 
case material O22662 Lbs/cw.iach 


The output for this example is provided in Table 


(IvV-2). 
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D. PROCEDURAL FLOWCHART 


Input data 


Calculate He eee drag of 

boost phase based on cruise 

speed drag and initial and 
final speeds 


Estimate thrust required for 
boost and calculate booster 
burn time 


Determine average altitude 
during boost and calculate 
ambient pressure 


Refine thrust required for 
boost and determine ace of 
booster propellan 


snc a 


Calculate pressure ratio based 
on optimum hamber ee 
to case weight relationship 


Calculate thrust coefficient, 
throat diameter and exit 
diameter 


Is exit diameter smaller than 
base diameter? 
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Lower chamber pressure 
Increase 2xit pressure 


serra anne 


Is (exit pressure/ambient 
. pressure) greater than 0.4? 


Is (exit pressure/a mb 
pressure) less than 


"e222 ewe @2e@ ee oe 3cr@e2 @e@e @e@e@ @e@e@ oo 


NP 


Set exit diameter equal to 
base diameter 


Lower chanber pressure 


Is the nozzl2= exit diameter 
larger than 95% of the base 
dianeter? 


Is the chamber Pecenure less 
than 100 PS 


Lower PLIM 
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Set) burn: rate equal: toy 1.25 
and calculate résulting web 
thixkness 


Determine burn area and 
required web thickness 


Are the two web thickness 
values egual? 


Calculate di nensions of the 
booster and nozzle 


Calculate required sustainer 
uro time 


Determine average altitude 
during cruis2= and calculate 
ambient pressure 


Calculate thrust required to 
offset cruise drag and 
calculate weight of sustainer 
propellant 


Determine mass flow rate and 
nozzle area ratio 


Calculate Mach of the exhaust 


Sa 
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Calculate nozzle pressure 
ratio and estimated thrust 
coefficient 


Calculate sustainer chamber 
pressure 


PST? 


eee less than 


NO 


Refine thrust coefficient 
Does thrust coefficient equal 
estimated thrust coefficient? 


Set them equal 


232 me @ @e@ @@@e@ e@@e Pee! Oe e@ @ @ 2 @ ee ew ow oe 


Is the burn area less than the 
motor cross sectional area? 


| Calculate web thickness 


-—fslenises sf. tiemest.___ 


Is the web thickness less than ~\yYeS 
86% of the motor radius? 


NO 


Lower th2 burn rate 


Set burn ar a equal to motor 
cross-section area and lower 
burn rate to correspond 


~~-#EED EELS Sp SEEESErEe 


Establish dimensions of sus- 
tainer motor 
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Determine average altitude 
during cruise and calculate 
ambient pressure 


Calculate thrust required to 
oftset crugse drag and 
calculate weight of slstainer 
propellant 


ewe e @ ee @2e ewe owe 2 @2ee 2 eee @w@e eo oe 


Calculate pressure ratio based 
on optimum hamber pressure 
to case weight relationship 


Calculate thrusti.coefficisant, 
throat diameter and exit 
diameter 


s exit diameter smaller than 
base diameter? 


Is chamber pressure larger 
than PLIM? 


Lower chamber pressure 
Increase 2xit pressure 


Is (exit pressure/ambient 
pressure) greater than 0.4? 


a ke 
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Is (exit pressure/ambient 
: pressure) Less than 2.5? 


NO 


Set 2xit diameter equal to 
base diameter 


sure less than 
LIM? 


Lower chamber pressure 


Is nozzle exit diameter larger 
than 95% of the base diameter? 


Is the chamber pressure less 
than.250..PS1? 


Is burn area smaller than 
motor cross-section area? 


Calculate required web thick- 
ness and web thickness 
resulting fron the burn rate 


pce oe ee_steh £2es. 
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Set burn area equal to motor 
cross-section area and 
calculate corresponding burn 
rate 


Calculate dimensions of sus- 


tainer motor and nozzle Gay 
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Ve. AERODYNAMIC COEFFICIENTS 


A. DESCRIPTION AND ORIGIN 

This program is the current edition of a program which 
originated at the Naval Ship Research and Development Center 
in 1971 (Ref. 5). The program was written in FORTRAN II for 
use on the IBM 7090 digital computer. It was developed as a 
method for predicting the static, longitudenal aerodynamic 
characteristics of typical missile configurations with the 
ccntrcl surfaces in a "plus" attitude. The original program 
computed the aerodynamic characteristics for missiles oper- 
ating at angles of attack up to 180 degrees. The effects of 
control surface deflections for all modes of aerodynamic 
ccntrcl are taken into account. The method was based on the 
then well known linear, nonlinear crossflow and slender kody 
theories with empirical modificaticas to provide the high 
angle of attack capability. 

The program was modified and, presumably, updated in 
1S74 by F. A. Kuster, Jr., of the Naval Air Development Cen- 
ter. In 1980, the program was modified for use on the Naval 
Postgraduate School IBM 360 computer system by D. Redmon 
(Ref. 1J. The current version of the program was modified 
fcr use on the new Naval Postgraduate School IBM 370 compu- 
ter system. It has been expanded to provide graphical pre- 
Sentation of the output data. 

It must be emphasized at this point that the current 
program edition is not in a completely finished state. 
Scmewhere in the history of the program after its initial 
establishment on the IBM 360, errors were introduced during 
the modifications. At present, these errors do not prevent 
the use of the program and the output data is considered to 
be correct for trend-observance purposes. 
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Specifically, the program does not produce any drag-rise 
phenomena for either the wings or the tails when Cd is 
cbserved as a function of Mach number. Additionally, the 
decline of the drag coefficient above Mach 1.0 is not smooth 
er as prolonged as is found experimentally. It is very fro- 
bable that these two failings of the program are linked toa 
ccmmon error inserted accidentally in the process of tailor- 
ing the program for use on the IBM 370. In order to 
temporarily smcoth over the graphical discontinuities, expo- 
nential decay functions were inserted. They are clearly 
marked in the the program listing for removal when the pro- 
gram is corrected. 

The input to the program is composed of a detailed list- 
ing of the dimensions of the missile to be analyzed. The 
current version of the program will consider a missile which 
has four symmetrical wings and four symmetrical tails. The 
missile may be either canard or tail equipped and either 
wing or canard or tail controlled. The program assumes the 
centrol surface is the tail, however, the input data is 
“mislabelled" to produce the proper configuration. For ins- 
tance, if the missile is a wing control missile, the wing 
data is input as the tail and the tail data as the wing. 

For a canard controlled missile, the canard data is input as 
the tail. Figures (V¥-1) and (V-3) show two typical missile 
configurations and where the input parameters for the pro- 
gram are obtained. 


Be. USER INSTRUCTIONS 

If it is desired to abort the operation of this program 
prematurely, two methods are available. When the program is 
waiting for data entry, press [ENTER]. When the program is 
hot waiting for data entry but is processing, type "HX" and 
press [ENTER] . Both actions will return the terminal to CMS 
mode. 
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When the screen becomes full, or "MORE...." appears in 
the status area, clear the screen by pressing j{ALT]| and 
simultaneously. At several points in the program, 
the user will be directed to "CLEAR SCREEN AND INPUT 0". 
This is to provide proper positioning of the output on the 
screen for ease of reading. If any other symbol than 
"GFA should appear in the lower left of the screen, 
press and continue. 

1. Turn the terminal on with the red Switch. 

2. When the large "NPS" appears after about 30 seconds, 
press (RESET], then press [ENTER]. 

3. When "CP READ" appears in the status area on the 
lower right of the screen, type "L nnnnoP", where nnnn is 
your 4-digit user number, then press [ENTER]. 

4. You will now be asked for your password. Type it in 
{the characters will not appear on the screen), then press 
(entza}. 

5. Your personal file must contain a PROFILE EXEC rout- 
ine with the appropriate Fortran GLOBAL statement. If it 
does not, type "GLOBAL TXTLIB FORTMOD2 MOD2EEH", then press 
(enraa}. 

aenecype. “LINK. 20 xxxxPioi9:1) AS 192 RR",. where xxxx..is 
the 4-digit user number for the project file, then press 
(evra). 

7. You will now be asked for the project file password. 
Type it in (the characters will not appear on the screen), 
then press [ENTER]. 

8. Type “ACCESS 192 B" and press [ENTER]. 


9. Press |ALT] and [CLEAR] simultaneously to clear 


screen. 


10. Type “"LAERO1" and press [ENTER]. 


11. Input the following data as it is requested on the 
screen by typing the data and then pressing |ENTER]. Ensure 
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that the data is input as either decimal or integer as spe- 
cified. The terms in parenthesis below are the program 


variable names. 


payee the following as integer values unless otherwise 
noted The integers must be two digit integers (three=03). 


Gcntrol constant (ICSC) Ol=-Tail control 
; O2-Wing control 
03-Canard control 


Nose shape (INOSE) 01-Ellipsoid 
Q02-Ogive 
Q3-Cone 


Number of control deflections (IDT) Less than 11. . 
You will now be asked for the control deflections in 
degrees, decimal values. 


Number of Mach numbers (IM) Less than 25 _. 
You will now be asked for the Mach numbers, decimal 
values. Each Mach number will produce a separate 
table and plot of output data. 


Number of angles of attack Gael Less than 25 . 
YOu will now be asked for the angles of attack in 
degrees, decimal values. 


Number of configurations, (NBODY No restrictions 
Each configaration will restart the program. Only 
the last configuration will produce the written 
output. 

Wing planform (ISWPW) Q1-Not delta 

Q2-Delta 

Wing position (IAFBW) ee oe after wing. 

01-No body after wing 

Wing sweep constant (ISWEPW) 00-Delta planform or 


Unswept leading 
edge 8 
01-Swept leading edge 
Number of wings (NWING) 04 


Tail planform (ISWPT) Q1-Not delta 
02-Delta 
Tail position ({IAFBT) 0Q-No bod arter tail 
Q1-Body after tail 
Tail sweep constant (ISWEPT) 0O0-Delta planform or 


Unswept leading 
edge : 
Q01-Swept leading edge 
Number of tails (NTAIL) 04. 
Input the following values as decimal numbers: 
Wing tip-to-chord ratio (XLAMW) 
Wing leading edge sweep (CLAMW) Degrees 
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Wing span including body (BA) 
Wing root chord (CROOTW) 
Exposed wing area, 2 panels (SW) 


Wing mean ee chord (XMACW) 


Distance from nose 
Wing leading edge (XWING ) 


Wing thickness-to-chord ratio (TOVCH) 


Tail tip-to-chord ratio (XLAMT) 
Tail leading edge sweep (CLAMT) 
Tail span including body (BT) 
fail Toot chord (CROOTT) . 
Exposed tail area, 2 panels (ST) 


Tail mean geometric chord (XMACT) 


Distance from nose 
tail leading eae (XTAIL) 


Tail thickness-to-chord ratio (TOVCT) 


Altitude (HT) 

Body diameter (D) 

Missile length (XL) 

Nose length (XLNOSE) 

Distance from nose to CG (XCG) 
Reference area (AREA) 
Reference length (XREF) 
Engine code (ENGINE) 


Inlet code (ENLET) 


Boat tail angle (BETA) 

Base diameter (DBASE) 
Nozzle exit diameter (DJET) 
Boat tail length (XLABOD) 
Protuberance drag (CDPROT) 


EE eee wg with experimental 
values, Reynolds number (REFT) 
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Feet 
Feet 
Square feet 
Feet 


Feet 


Degrees 
Feet 

Feet 
Square feet 
Feet 


Feet 


Feet 
Feet 
Feet 
Feet 
Feet 
Square feet 
Feet 


Q0.0-Turbofan 
1.0-Rocket 


Q0.0-Flush 
1.0-Extended 


Degrees 

Feet 

Feet 

Feet 

(Coefficient value) 


(Dimensionless) 


12. If you desire to rerun the problem, or want to run a 
new problem, answer the questions appropriately when asked 
by the terminal after the execution of the current problem. 

13. To receive the printout and plot of your encounters, 
answer "no" to rerunning or restarting the problem when 
asked by the terminal and follow the directions presented on 
the screen. 

14. Upon completion of the program, type "LOGOFF" and 


Switch. 


15. Turn the terminal off with the red 


C. EXAMPLE PROBLEMS 
Table (V-5) identifies the output variables as they 
appear in the output tables. 
i Example Va. Tail control missile 
Figure (V-1) illustrates the missile used in this 
example. The dimensions for this missile and cther input 
parameters are contained in Table (¥-1). The output is 
Shown in Table (V-2) and Pigure (V-2). 
2. Example V-It. Canard control missile 
Figure (V-3) illustrates the conard configuration 
missile used in this example. The input data is contained 
in Table (V-3). The output is displayed in Table (V-4) and 
Figure (V-4). 


72 


Ss sos of ffdev <“G ,@t-(aes@ eo? BezSes oss? 
Loins oot) Yloteiseortege ent > izeug" eee 


aoe a ne 20 SW ris ) We? 


XLNOSE 


Figure (V¥-1). Tail control missile as used in Example V-A 


13 


sn wl@esed a) heey cx siiselw 
od > 


3 


, 


= 7 


THE FCLLChING TABLE CONTAINS THE INPUT DATA FOR 


TAIL CCNTROL MISSILE 


EXAMPLE V-A. 


DOONMNNIWDOODDOVOCTIOOOOrFG oo 
DVOOOFr-OMOOCOCOFINNFOOONONDG Own 
CODDMDORTOVOP sr TONFOCHORNHFHTOOCOLTNOO 
SNH Ota H4 HS Oe eee eee to reer esve te te eevee er eee ee 
- OF ANNAHOON AODDONDOCOVOMONDHDON00090 
~ + (2) 
>) 
ee ee oO 
Ww uw ina) 
oO [e) 
a io) 
IU uu 
ee 
Q YO oO 
aw oe z Zz 
dee (a) — oe — 
Zw ae GY it eS) 
<i —_ mm 
oMm =e WwW oe aU) 
to =) e J _- 
laa] a oe ve ry 
faa) Wee (9 Alameda oe O ~ — 
= =Oo- woe _- (a) ot tw 
Oe eo WU eee HUY Ww & = U. mee 
aw OdaFr Ui z Ww WY ee w eeve (va) Oa 
“> ZreYM ett w 2 WO « Wo (eo) ow 
Zeaee —> Oil <a 4 il ce Dee J t) oo a) Fd as 
won =SZCW aS Reus 8 De Ree - ke HZ 
NOZ on jaye) ous MeeIetID Nee tt oe ow 
WZ ev veegDm - aQO ef 29OQ ~YyOZanw~ dIOnMO oO eh oo 
ea Ye WU OL i OO Ot WI «Tt aw ax<— 
SA Fe UZ O JN BAe DZrORPSeaewe ZFOrFOD uw wn 
—eD <OW ~O WH~A w>nstr ¢ wW> AT eX elu 
AWW eR Oer ww OF W LWOQ AVUOOUWO AVLDOO Cs zou 
me DJ ee be 1 Or Jb OzOr [5 ol Ge Na-K ae ] -r Oo —q oy We 
WLUNAAIANZEO er) Mees (co nd oo © — US Sey Oe ™ ee 
ieee’ lc Aas Ak © Ade © ~rew OO cw — O UW seed 
JOWLD Fee Cl ol lo] mld FOND WYO ot YNO Ba | moe +I 
eb DMO ADL NO 9ZOMR DEE O-7B~R DE W °rT~ uUt-ao 
rw SZAOZqANGYZANADH SL we <A we (eo) a Fe 1%) bee ao} 
ZU IDAL SW DO -ZWO UOT ON ULL IAI YN eT ed ee «© 
Ive 2W2Z FROZE R OY DIWYFWN DIWOZW ee TOON ee et 


Ee —feah 4b" a OL 9 3 St ts Di Da St 
NeeZrliyoaaapr Bil SU OS Fo P= Char) C= Ge Na rrecilg Gos OC Ut Vin ied 
2.0 22 ke ae mH (LO (ete IL Zeze u< 
OOO044UI>WWL I> WWE EI4COz UTO~WO4A2Z vo 
Ot FT OQWOFPROWVOFBOYOA DEAUHOLt Owego 


Os alas IU Use eh ea NIT tates HI Hees Yom KU J 


EDREOctii en zr = 
= Dit S991 Ue Sit 
x'9Z0 e Leese 


WZ LU WW ee Det x TU 


Te u&u QOMW 


el ACU CI eee oe eon = 
YU Oo we OGaw]a2 4aDue- ao 


agonn 


PwWeyewOTIMerdawOoT FWY UOW + Ur-O 
Sead Ca ren CU Se RU een <1) (es) 


Ro a = GS) QA We WOW rR We 


EY Se OS) SAVED ti UE DiGs Sis Math Say a Se ee ee 
FNWYWETSHOR FS TAOFAZMNSW ZN2L Naa NH SO TI Z DEN 


ZUMwMIYN « 
UZAzADOSFZHa0N 


ON 
MU JO TIL £ 
OfA02zaDnNAaIY 
tt at ts od et 
ee ee ewe eer ee 


ee ee See 


ANOIFNVE OO 


me | 
<HoO 
_ 
ob 
Orra 
Zzadmu 
— Fu, Fk 
nan 
Zee 
et eee ee 


= =.= 
QANMS 
or) 


=> le Lee IT et ea FR La zat aoo 
Ze tO O FOZRRY OFDR eS YP UNWIU HUMmaD 

mmm Ee at Ke OO We = OF 
Af FSW CY ERR Rr WR RIN Re ol er 
oP a 9 VYVZuz> O OF OO DBD O-aAwewaneoe 
A494 de— OK te eottHKOrdaat or HK ZOU et 
RI £Y CF EVI ISR oe ee OF mK str ee DIVYHO 
SMUNDOUN XR MUON XK KR LOX KK ef MW MOD KU 


——— ee eee 


re nn me me eee eee eee ree One ese Oe eee eS 
TNDO™ DROANM FNOMDIOANYDTIMN OM DROANMGS 
ot St SIN OI NIN AION NAIC MOOT EA INCI FS TPT 


DELTA 


ALPHA 


Input data for Example V-A 


Table (V-1). 


74 


CLbW CLTP CLBT CLB COI CNWP CNIP CLTO COTD CN CA XCPW XCPT XCPB KCP Cm 


CLWP 


COV TNO 
OTe TNR ANS 


NAPOTONSO 
' 


tt |e 


Paton’ 


NW Dae AAO ao 
SF Oornwuee 


eoereeeeeve 


a laalaaloatan Loalon Dnalaa] 


EMEC AMM 
me=MAMSaOCO 


esreee rete 


Ont et at ee NY 


et et et et st tee 
| el aed eel peel el ed od 
oeeeeee ee 


DWM ALA ALA AS 


lelelelololelololqo} 
TRIO ATFAROHK 
eececeeveeees 


AMMO AA een 


AOCRDAS ME 
SIADDODSCAHE 


eeoevreeveeee 


ococoessco 


NOONKDA Ss 
WF ANNCIMSA 
esve@seereeeee 


ANON THAOE ON 


coooo0ececd 
cooocoose 
ewer resnee 


ooooooosse 


lofololololololajc} 
>WVOCOCC S00 
eecverteeene 


coocooceoses 


INTRONS EN 
Coonan tnt 
eecverereeen 


onoococooo 
' 


’ 


DOWMNDAARASY 
NPANDANAE 
erveveeee 


AHN AMGSFAMO 


PFOOMaMAFNOO 
OA-NAMAPE OMENS 
ee 


eee ee ON 


oc C0O9090 =m > 


ee 


oo 


OM OLHABONOH tt 
COS TDHSNNAL CH AOe 


ose e eee 2 20D 


oo ocpCcocoeeSdca + 


vo 
a 


OMMASDOMNAGAONW VU 
99 CO Set NN 


cee eee eee N 


oo cocaescocoon 


~ 


coo 


OM FNUMOFONOCO! °& 
CO COO RH m rN rO 


ee 


alo tololololalolololo]. 45) 
an 


oe 


SOF TOPONALTOF 
Da aCIMPNOE DH O 


ceereseereee WY 


00 900990099500 


a 
oa 


On FONBOLPMANN os 
DDO NDNHDANTOO— 


O09 HH NMMOS ANO-O 


ww 
a od 


NA EANONORMONOO 
FR PO—NPOYNOOG 


eevee ee ee WY 


30 Ona mts 


OF INANR SSK NMOT OM . 


eevee er ede N 


O90 HNNPTNOE DAU O " 


eoeee 
Sa NAST AOE DACOY 


a" 
== 


© oma 
“OY 


Table (V-2). 


a0 

N i-i=) 
r oO 

OM NVNDODNsHOD °o 


M 


XCPT = XCPB xCP 


CA XCPW 


CN 


CLTO COTO 


CLbT CLB COI CNWP CNTP 


CLBW CLTP 


CLwP 


Tcl COTOT 


_ 


NETO MNNTOR OO 
ONO AK TONY NHNOD 
COG) Oia: 7O\@) 168) 8) 1S) 
DODO Daw 
DO et et 

CSU est iL 90 


oFmsmramoronst 
eM belo heed ae Ler =) 
CG. O20) ' a) .8)'S 26' 6) 1s 8 


ALT PP perry 


COANE NOMA M = 
OWCHe DAE SOOC 
ceereeeteeee & 


loleleor re tt ia ia) 


et es tt ee ot ee 
ld ial td a ld 
O88! 0) 6.6) (0.6, .@, 0. :@ 


Walabalral al ar al alata) a) 


elolelolelajojclalolo} 
AAG AAHAAMHH 
Ce Gee oS +e) 8) 610 0. @ 


AAMAMAMA MMMM MN 


TAMDADOR ANE 
FROME RHROOs 


jolelelolelololor lar.) 


PNONAARANA HA 
SAMA OO aN 
@) 68) @ B19 0 8 0 8:0 


SHaNATNVEOTSO 
~ 


dt et tt et 
lelolololelelelolejlojlo} 
ceeoeesereerere 


wlolelelalolololololo) 


AAGSARAVAOOOO 
Wolval'alvalval ar al al al'al a) 
oe eet eoeee ere 


lolejelelolololololel—=) 


oODrNnamaooOFo ra 
DNNSOOE Ot cM 
eeeeteoer eves 


Lolo elelal@lolal lor) 


ONNONNOADA 
DIANANANOMNSS 
cee ereeeeeee 


DOM RKN ASFA 


APOPKHE OM FONDS 
DOM =ALORNIQONG 
eer eee eee Om 
lolelelelelololor Tal. bales 

ee 


oo 
SOMOPHAOON OH It 
CODTAaNNMNOS Ome 
ceore eee eo e CO 
lolol elololelojslalololor 

oo 

io} 
FaPFORALPOD—M O 
ANNNNM NMS et 
eeeteeeeer ee N 
lolol elalelolelolelojlojlo)'s) 

“~ 

oo 
CDOs KONNOEUE I « 
DNNDOE DBAONM&HO 
eee ete eve eve 2 en 


lololelelolelojer Lala 45) 
an 


On 
SOTA DOFONALTOLE 
OSM NOE DD 2 
seer eeeve eee YW 


oo9aqc000%7o0S9 


om 
CHA TONAOLTANN 9 
CIDNODNAHDANACO 
eeeeevre eee eo eit 
ON SOUND Ie 

ui 

an 
TENT oOeOrsnosd 
VOOaeN SF ODVeNPOO 
oor evee eee ee iy 
DO me HH ON 


vr 
FIPADVNONDO~-OO 
EFAMaSUOOU0O—N em 
eeeeeeeree es oO *@ 
VSEINMTNOE DEO NO 
ata 


“ 


XCPT xCP6 XCP c 


XCPW 


CA 


CN 


CLTO COTO 


CETP ELAT CLB COI CNWP CNTP 


CLbwW 


7 
EFOCASNOSNMOM 
MNSFOSHATIO 
Sid: 9!) ,@! © 8) @- 0:6) ©. 6 
CaNANAE DONW.E 
me ey NNO 
ESCs cheat ENC) 


DHE NOOWNDO 
mae eS NN 
eoeertevee eee 


ANF PEF PPPS 


Cn AE NOMMM 
ONWDemMwes™ DAO 
oere eet et eon 


COO we wt HUN 


te yt ty ed et ed 
ld dd ee ld 
eeteeeeeeee 
Vala al alalalal al alata) 


loleolelalolalolelelelo) 
PRAPARDAGANH 
eereeerese ee 


AMAAMAMAMMA oy 


CANONS AOMOLP 
CONDDSCCO SNS 
eoeroeeeeeveeoe 
te Ot ed ed ed ond ed ell 


FOP LONNANAIA 
WWABORTFOTIAN 
eoeeee0nreeee 
ANAS EN OF Con 

— 


AAA MAA ACM Mc cy 
lelelelejleloleloleleya) 
eeereevever ee @ 
elolsloleolalolelolola) 


Oe CONnNWSES 
et et et ed at ed ed ed 
eer ee eevee 
eh tt ed ad ed et ed ed 


NIN ON SH SMO 
NNNNAMSNOPDO 
oeeeeeeeeee 
et tt tt tt 


SCNNOANOAAD— 
SINDAANDOANA RS 
eoecereteoeeseve 


SOK AHNAMANTIANS 


DODAR IFAM ARALOO 
Qa VN POTH TOMAS 
ee ee tee eee omy 
lolelelelelelor TTI. haley) 
ee 

so 
OMODTaHARONO 1 
QOOCKNUNMDNOFAae 
eeee eee eee ecm 
lalelelalslololclololap 4 
[ST] 

jon) 


PONMEANSORS OO 
MFT PIOMNYNO 
oe ee eeereeree 
cocoococoocooon 
Lied 
oo 
NOW -OONTONnH @ 
NNNNATOOEDOrO 
ore er eee ee & oO 
a th et ng mt OO 
an 
Q~ 
COST DOrONALTIOL 
AMINA TAO-= OP GC 
eoeereereeeeere Y 
OADQQNCF09009D090 o 


on 
Ca FOMAOLMUIN Oy 
PFIONDIAKDANA DOOR 
ceeceoe eee eee oy ® 
QDS ANAGSADHRO 
wa 
Ie 
OWMAPFHADNOSL-FROQO 
OVUHKINAWO=TFONUg 
oeeevereeoree Wy 
tt tt INO 
N 


vv 
APFNOAMTOATOOO 
NAO DWODP Pe oy 
coeoee eevee es 0D @ 
KNOTT NOOO 

“~~ 

<t 

ee ee 8 8 6 8 8 © Ome 
QDaINAFAVECADIAGQ 
“OY 


Jutput of Example V-A 


Ths) 


“ Tey 


CA XCPW XCPT XCPB xCP CM 


CN 


CtLTO COTD 


CNWP CNTP 


Col 


CLB 


VN NOTORPTS 
NFOT OM HHO RN 
oe ee eee esos 
KWON FONSTCDON 
SONNE 
(beth Ue iGe De et Ge Tha} 


DOOEMAMMOALTO-~ 
MOIS OAT AMIN 
OO) SF OAS (6:18) @ 8-0 


Se i a 


ORE AENAC AMM 
CWDaMAK CPO 
a8) 6) (9 (6) OF 6 '@' 2 0 ie 


COOK aa INN 


et ss tee 
~~~ RR eRe 
eee ree ere ee 


Palalal-al al al al al al ala) 


cooecocoocea 
PAPA AOAAA HD 
*ee@ewreeeteeee 
AMAMAMAMAAAA NES 


BONWNEONSCUNMN 
CONS RINN AO 
o*#e eee eeeee 
me ts so 


QNOKHKaAmMomuS 
Na PDOerrer DEO 
oeeeeee eee @ 
NAAT NO~ OPON 

——_ 


WODODODODNOOONW 
elolololololalolalolo) 
eee ®@eeeeee 


ocoooco0s00090 


PONS tOOOr~ OW 
Ll ed Dl ah ed ed? OOD ®) 
eeeeveeeteee@ 
et tt OS ed ad ed ng ed 


MANS ROMA OME 
DBOMortoHnmnd 
OPO) ORI CY SOI) 
at met III CIOIOUN 


. 


SBNNOCWNOSA A 
OCOCNAANTANDS 
eeoereeeee eee 


DOH SNAG FPOBO 


SwDBnNCHDeODrFvO 
ASIN SNOT ME oe 
eeete eee oe 8 oe OMAN 
S99 SOON aa NNO 

oe 


oo 
DAD LPH ARONO— ii # 
SCO HKHNNMDNOFO Oe 
eerteeerereev es 20D 
elolelololojolololololar 4 
vo 
j=] 
NDDONSTEAHMA~ UO 
WNOOO DOE Eee 
(eters so Se" se" y 8 N 
lelolelelelelelolejlelejlo}'s) 
_ 
co 
ANFOMNMTINDNO HW * 
DBBDODVDTO$NMNOKO 
oeF ec eee ee © oil 
ttt SI HIE NINN OO 
an 


c= 
BOSS ODE 0MNMPoOT 
Oma NMPNOF DR 2 


eoeotreereeeeee YY 
lofololol als lalelalola) m 
on 


DHF DMPOLANG os 
SOODNDA-DAND OOS 
eeeeeersteeeele 
CORN SANOFO 

wit 


te 
DO ON~MOF FONDS 
Va FROME =WOO 
eeerer eee eee VY 
tt ot St INNA 
N 
vr 
OTN ONE NOOO 
Ae AO™-OUO— OO em 
eeeteevee eee 0D @ 
NANT NOP ORONO 
~—mt 
= s 
2 t=O 
FosSo 
=U 


O. 
l. 
2. 
3. 
4e 
Ss 
6. 
Up 
ue 


Table (V-2). 


76 


OFTTNCRHOO-TWNO 


FUME OWNAKIOADRO 
re) Oo el -@ Leyes 8 Oe e 08 
O-aNISIA™ BROAISEM 
NNN NNNN OCHS 
(Pe SPS WaT SOU WSL Cees a] 
a DDO OR ONOO 
Oo AAT OUNe TD 
x< essere e eee ee 
DIAS FP PPP Pe 
OD DeEorNnOMMAM— 
a ONWTeMEDOOs 
(jo) seer everer ee ee 
K COCO N RH eS HNN 
al et tt et St et ot 
[oN Nd ad ad ad od al 
(=) ove (e" 27 6 8:10" ee Be 
xX LAMINA DIAIAIA WIAA 
Ez DOCCCOCCOCSCSoSO 
A APAPSSAAAAS 
Ww eereeee eee @ 
K AMANSAMMAAMMMA 
I CMOSNENDOAW 
VO ANNNNGM AST ONOPE 
oeee eevee eee 
td at ed ed et ed ee 
ZF ANDNSTOARNOEaA 
(eel alliatd'al at cla’ 
seer eee er oee 
MAMFNOE DROWN 
ot mt 
CG BSOOCOCOOAAAGN 
— RH Ht HOC OO 
S orev eeevesvee 
Oo 9000000900090 
0D OWT NRK OOENM 
KF NNNNNN a ete 
23 oer ere ee eee 
YQ NNNNNNNNNN 
a2 AMNOFOCODOME 
Ee SPT ONSCE DON 
za oesveeeee ee ve 
VY NNNNNANNNMGOA SS 
Qa OCONNOCWNNOOSA« 
ZF DONCANTANORK 
z oeenereveeeee ee 
QD CDOK—NAMRSNWGO 
= SD ~O0DDNFFMDAVO 
QD ant n=Onvoranm 
oO eee eee eee & OMA 
SASSI SAAN NTA 
ee 
oo 
2 OMDMTAaTOONOAAH 
a4 TOCANN NNO FOr 
Yu oeerevreereeose OM 
lelelolololojojlelolololas 4 
Va 
a 
= NPOTDONFODON VU 
2 RFK DODD ORAD 
— oeevt®eeeveevree 
Y DOVCOCOOCOCOCOCON 
~ 
oo 
QAP eOrromnon * 
— FPP POOF Oe HNrKO 
med eee tse eo ee eOu 
QO NNNNNNANNNS SO 
aw 
a~m 
ZF DODACDOKONALTOCE 
CONN POOF OF AQ 
a eeevereeeeoeee WD 
Y annosoo0ctessS 
ona 
Qi CH TOMAOITMANN *> 
FZ OCONDANSHDNHVUA 0O-4 
S| eereeree er ee cif _ @ 
VY COAaANMATON ORO 
un 
in 
= DSONTIDNANNO—O0 
D0 NATANFOMWVONOVG 
-_ eesveee ee eee Vy 
990 AM6 SRINNINIOS 
alors) ™ 
oe rT 
QO0F NHKVNDUMaM=TTOO 
cat DEAT IP TOGO OM 
Nn eoeoeeevree ese? 20 0 
fa MNT NOE DYVO=N HO 
=e <= 
Cad o08e© © @© 8 © 6 8 8 Oem OH) 
dv OR-NATN OF DOFocNd 
maa “OU 
(Continued) 


ipoutraday 
m1 


- 
7 


SH at pag 


! 


1 
' 
BS See A Stitt | SSS SSS 4 


i 
' 
r 
' 
i 
+ 
i 
i 
' 
a 
i 
i 
' 
if 
i 
i 
! 
! 


EO a 


lee tarts ieee 
' 
' 
‘ 
' 
espe 


Y Pao seer TA ae ee ee lS ee Teroeeere a i (Se c fe 

' ' i i i i i i 

' i i i i i i i i 

i i i i i i i i - 

i i i i a i i i *° 

mee cc ha © aaa caer a a) ae SK Sees Gas Gia oe en ee 

\ i i i é i i i 

' i i i ' ' i ' i 

' ! ' ! ' i] i] ' 

i i i i ' { i i g 

mao = ie feo amo ara NR Ee [a oa as, 

i i \ ] i 

' ' ' 1 1 ! 

i ' i i i 3 

i i i t i 8 

(Fg ean a opi pe Sea, ake sence | Sea a aaa ee, le a Cee eC) 

i] ‘ 1 ! ! 

i i i ' 

' ! 1 ! 

i i i i 1 8 

[a SI ea decent | a Lem ly a  etenies (fig miami weal ated Gimin. <tes <a ee, © 

i i i i i 

i i i i | 

! i] ' ! t ' 

i i i i i 8 
eo’o eo-af- G ¢0'0! oo ) a9 oo of 

w ul Hs 


Output data plot for Example V-A 
77 


W=2)). 


Figure 


h-¥ =tguaes ‘se? w6hg ands Se 


- a 


XCG 


. 


Figure (V-3). 


XLNOSE 


D,XREF 


| 

| 
| CROOTW 
| 


BETA ae SAEs 


ores 


Canard control missile for Example V-B 


78 


THE FCLLOWING TABLE CONTAINS THE INPUT OATA FUR 


CANARD CCNTROL MISSILE 


EXAMPLE V-6. 


MVD AHHH 4 
~— 


Ww 
Le) 
(a) 
w 
ee 
Q D) 
oa Oars 
ave C3) "ea 
zw Fa es) 
cia —~ <a 
Sa) me ey 
tO 4 
i (a4 
fon] War 
. eon 
Oe ook ut 
zw Oars 
~> Zz eM 
“Frsaee > 
now rTOW4 
NOZ (ate) 
wo oor ~ 


aww FeHOr w 
ROI err HOR 
NQuUMtauzZ 9) 
tale Ul o Sutera Sa I Goad 
AQWLD Free 


se Y MODODD (ND ND tm DOSE OEE OO EE OK 
RWIS MAY ZIGGY SZANTIOMY We ater ws 
ZntCDdNth— ITU AT StU Yoo s AS wc as 4 


Pons —_ 


LEADING cOGE: 


TOCOONOGCOACOFIMORVOCDOMO el alele) 
HAOOTENOWMTOODAOFODROOCOOOOND oo000 
MUN DADQWNOTWOCAODINGEMOCOKCNOWNDIOCLIONO 


wre eee eee te eteeeee tet eee eoevee eevee 
OE NSFONDTONNODOONCOONNEFOORHONOOHO 
iva) faa) oOo - 
oy 
© 
~~ 
ee 
eas ve ae 
* O > 
_~ (Va) ee tu 
YY Fe Ww uJ NZoe 
wy i@9 ese UL eee WY ON 
Ue a we: a wo oO Ouw 
Cees! Iu 0D oe J hne Ps Yo 
3) De ter -~ - Wa 
Ow YAeOAtyuO Yew Jt 2 OW 
ZO >+AAZXee~ AO~e¢ (>) ob— ee 
tates (QW ae. eT UY eyed mx aot 
BQA2?O2ZC2DAe ZFOrD uw Wn 
USO ne W> <r oi Ud 
LWOH AYNOOUWOQ ALOO = zow 
ok 70) —-LOzO -r oO <I oC oe 
YAsstoo wo WM Oo “ war ° 
(‘D) ~taw © Zrew -e Se ise 
SUD FASANO WY) oop YHOO a4 amem -+-IirI<_ 
QO wexytT~— Wey 
2) 57%) eA ke) 


Seek —DUie TZ 


daw Z7wWZEROIZEE Ie DDworw SIwWornwn ee TON Ieee 


(—S Ae ies ue 


mm UL 


Cat ED Nice eee FAS i AAU Nien Pa Via Yeo Ye Jie Je AT Hee HH] dente OU Sls Susie) 


Wee ZAG) Ae Oo ee ih eee OO) SA Or ie =O) p TOs pa | A Nn ANY I Sl a z 


we) es 16. 


lod Ja 


J34.272R SS Ae de Id ww we) Stacod Ye Lie 


OOOTTULISWWLIKWW TOMO SZ CLTO-O4IZ oO CO2ZzdD O; Be heal lo < 
Ot TT AWOTPFQDWOATEDVA OF ATW SOT Ont WWE UUern Ut<etw 


Ge vu ADAM 


AINA 


iw eet Du emery T RWI WOW eer oO 


Yo eee eas ea) ples ate al | all ea ae) id Se Soe nat of aT oa If ae mel I ea) Tp Neel eC cde i) 
O OO) wie paws aD <xco Cot he EGS) QA WH WOW RRR WR 
NY Si DVD Marae tice Valea O29) Sea a Sy et ea ere We SIO 
RPAUYNWTZAZSHORFSAIDITA ZTNSWZNFA HN AH DN OI EON OY 


ZOMOWO DWH se 
O7FZADEZAO 
By SxS) = 
—Z25D2---= 
wl ££ x» £ FO 
Onn =—JAa oW 
NOK H-SISL = 
YVYVZQ2ADBDMAMN 
et ty a es Pt 


me ee et es ee ee 


ANH TNOM OD 


Du we 
z-0 

od 
OFrFa 
zamw 
“Fi = 
ZMaIM 
Zee 


Sy 
CANM 
Cente kamen Rae 


Je SY aes eT ND ae OR LEZ IMINO 
Z— FY DED Y DeEORKN WN ww “INDY 
deotee Had EK Nae MYO ND Zam 
mm OK aT Kee OO uum —~ OE 
Ba Di I I SSO A ll ad col OS] eed | eo OS Df el Pl eal OU ee] OS Lee) 
=—2> OF OFUZTZE YO VO DM VARALewayrog 
Static Ona S at OKA eee RH SSW Ie wet 
eds fe £2 SOY MS Or em IU I CZ ZWwDIYO 
ZK IUDION KK KUDOM™ KE LOR KM AKL DID KO 


ee ee ee ee ee 


SN nnn ee ee Oe ee ee een ean ae eee 
TWOP DPIOANMISTONOH DOOAHINYDTN OF VDHROANMT 
St tt IOI ONION IOI IN IN 0999 OF OED TVET OV OD SE SP SP 


MACH 


DELTA 


ALPHA 


Input data for Example V-B 


Table (V-3). 


79 


CM 


xCP 


XCPT xXCPB 


CA XCPW 


CLYO coTOo CN 


CDI CNWP CNTP 


CLAW CLTP CLBT CLB 


CL WP 


SF CoromyunsNna 


o~ > 


toraerromn 
. 
1 


LE aannnan-o 
Tormrennanrg 


ee 
Lhd 


eaevrteeeve 
dh a i Dall el 


ON oman connn 
CO Rar athony 


OO OOOO OTOEC 


OF NaN mms 


SS coecoosuco 
NNNAINNNNG 


overt e et eve 
MAMAN AMG MSM 


Dalat ateatealoatealaaloaloal 
et et et et 


—— 
ee 


SH gpa gnaacane 


DD nor amnaao 


ry 
ee 
lola) 


rrr 
eeorree eee 


ooco°o0oco 


S28 nmvonnomsro 
ROONDOMATS 


° 
9° BaNnNim Ano 


22 gaoc 000000 


oo 
ee 


elelololelalojeleo) 


29 go0000000 


2° ga0000000 


oc 


fol Sal Sel elojola} 


SS gc0090000 


lelolololol@ialolol. 45) 


OW OSE NOM 
SANNA Sw 
eevee 


ooo700000 


SOSNENTADO 
wOoOLtmaman 
2S) '@) 10040) 8; Sie 
Oma INNAMSST 


NNO SD PMOAN 
DOnRKNMNAE ODS 
eer ee ee ear 
o0q000T90900 s+ 


ee 
oo 


APUOAMAOr Ph WH 
CaANMeAnn~aocar 
P50 S' £6) B48 8 


oooo0co°00o—o+£ 


eeI0 


Vo 
=] 


anny rnono™~ Oo 
loleloleloloiolola) 
eee ev ees ee 


alolelelelelalol olor.) 


wn 


oo 
SMNONOCD HEA « 
SaNAU HAST THEO 


an 
om 


NMAMOEANSS rot 
HUMANA TNO 3 
eoereevreere ee WY 


enoo0n7o0osr 


a) 
ON 


DOVer—O—M Om 
COL TNOMIANO— 
ae) @ (Ba) :0) je 81 aR oO. 
CMS CIN AS Pee Oo 


wit 
te 


=ANPFODODNNHNOCSD 
WANOSDONTOD 
erveee eee J 


90°0°°9°0--—— 


cond so 


Dal oe 
NAC O~DVNONM 
DIANITANGOE om . 
eeee ee ee 0) © ey 


AANNMPWOUVO HD " 


eee ee © © © ome Ss 


NALTOS™OPCOO 


su 
ze: 


“OU 


Table (V-4). 


XCPB XCP CM 


CA XCPW XCPT 


CN 


CLTO COTO 


COI CNWP CNTP 


CLB 


HNO KFOMNAOMS 
CDN INC MADMHHO4 
Dy Oy GO; Oe (2.2) 8) 2. 9 
Po alelelolalal tl. 

Ue Cit aS | 


MAN D-- PHONO 
Noe rn0oRere 
ode @¢ o 8 Os 6 
NOPE ERR REE 


ONOMOEAORNIN 
COR Mr =TrPFone 
ee) 0.8) 0 € | le © 16 @ 


CHKKNNMMON ITT 


C09009000000 
AINNIN NAN 
® @.8 8) 'e ©. @ 8 ©. @ «6 
AMM MAA MM MMe 


AAA AAA AM 
aa es se ee 
eeoereeeeveee 


APADAOAPACAAMA 


OADOL$DBDAVNCH 
wnrrsremasanmN 
O20: O80) @ 18;.0 oS 


lolololelslolojolololo) 


aon nN 
NE NET —oMn=0 
"8:19 19 PB 1626 BB Bs 


CO SK SNAATHNOF 


coocooo0o0oooco 
cooocococscosoo 


a’ 0) 2's: 0; 6 6 0.80 2 


ocoooso0ocseo 


NANNING 
tl et et tt 
eoeeteeeetoees 


lolelalelolololojlclola) 


ND PFONKEMOLOC 
SAHNMA SPN NOPE 
® ©. 2.6 6 8 € @ e's, .0 


eaoqcosecoecoec 


COA er aoNeKnrnE 
OmMoorrmonnnm 
eseeeeeeve eer 


00 Om SS NMOS 


DAMON ONSRNAN 
PFO DOs NINE ODL 
ee eee 6 © © 8 OPH 
ielolololalpigiolole€ lak J 

ee 


5a 
OF ANMAMN MO mmr it it 
COSA TAyDCoEe 
eeeeee eee © TQ 
lelelololelolelelolo). lar. 
vo 
a 
EEDAQCHNMMS VO 
loleleler tt lt) 
eorfe ee eee eee 
lolelelolelolololalalelals | 
on 
oo 
ANODLARPTOVAIEMA © 
at NMS TN OD RKO 
ee ee © eee ee oH 


elelolmlolelololeloloh ae) 
an 


Que 
SCBNMOMPMNNAELOT 
ACANANSAWMOS 3 
ep eeeesvereree VY 
ASBaDI7Oo590e0a0 


On 
ONBDOOAPLHOON et 
27 0OrD NO DNS 
ee eereetoeoevreeit @ 
SAC eA eAINGE ATO 

wn 

Tt 
SON DH TDR HOO 
MMWANO DO FONMOUO 
eeeereree eee e LY 
990790900900°7°-=4— 


= 


_ 
SOLlrwroecnvnownom 
NEN OT HENNA OM 
oe t@eeoee eee erme 
OO AHN SIN DO it O 
Jil 
ze 
a 
GaIArA 
aU 


U. 
1. 
2e 
3. 
te 
5. 
6. 
Te 
8. 


CA XCPW XCPT 


CN 


CLIDdD COTO 


CLIP CLBT CLB COl CNWP CNTP 


CLAW 


COOM sNrK— reco 
“OOM DOMO- 
Sie rie. 0l Ble: Sie Vee 8 
SS et ot et HOC OM 

' 


NANTONMOANAN 
WOOO UAT AML 
eeeeteeeeee 
AnNoDDOrRee~~e 


CNH OMAR-DOKNN 
COrRMR rons 
a0. 0) [0 Ci @ 9 @ 6,8) © 


CHS NNMMOATl SL 


cacootcoceo 
NIN QLIINN NOC 
Cy CR CCC rat Per rt 


AAMMOMM AMAA 


SOMA RA MMM oy 
et td tt ot 
eee eee eeseoe 


APARAAAAGADSS 


SOPRANSE INO 
WA Pt rr mma 
eoereeseosene 


wlolelelolelolelelolo) 


ANEOCOOsoC Oo 
FRLOSNOOM #N 
eeereeeeeese 


SOKNNGA ST SCHOO 


soo000090sc000 
ococ700scceo 
oeeteeee eee e 


cooo0caqc09nscos 


rrr saan 
NN ANRA NNN 


olololo) Sleleololelaro) 


lotee teal ot ln Baelenl al lad 
NAMSANDOO™ DO 
eeete ete eeeree 


aqqocos%ocoo 


CONSE Dae eAnTrE 
onvoranmananm 
BB, S.C @ VO. @: 0'0, 0 
C00 NINE 


ON PSDYINNANMOE AN 
SCOOQOmMNNMSPODOD SF 
eoeoee eee eee orth 
wlaisleleleleleloliel lend 
Ao 
SOLAMAMAMOMEK I 
BOCaANATASOSCDe 
eee tee ees 0D 
leolololelololelololor [ae 4 
Vo 

ES 


FLNOOFDOASCO vo 
he ht eg et IOI 
oe ete te eevee & 
aqooo0sec00seccon 
Lal 
co 
MARPLOONDLOO!N * 
NAMATANOO- DDO 
eeoee ete eee e ett 
Ssococoo0030900Cr7U 
an 
a= 
SOMNMOPENNUAE trot 
DIAN AMPANS™ A 
eevee reeeeere WY 
AAOSOAIAAATOS 
uw 
AN 
MONOD LADO om 
DOPLPDAGTNOM Os 
@ereeeese ey @ 
CPO ARS aA THD 
wi 
<—— 
OA PODMHKHNANNS 3OC 
DUIANAD OF DOKMNOD 
sere eet eevee mS) 
39099990—4— 
-_ 
aa 
ASN FONNANTOON 
TATOO ONSEUM os 
eeteeteeeee oDe 


SOHKHNNM TST OOF IO 


ee*e ee 8 © & © ome 


OeNATNOEDOODO 
=~OVU 


Jutput of Example V-B 


80 


*F wloasas To SUGIRE ) joy 


. 


’ 


— 


CA XCPW XCPT XCPB xce CM 


CN 


CLTID COTD 


CL.TP CLBT CLB CDI CNWP CNTP 


LW 


c 
7 


Man OL NON AR 
Neate neT eo AmMO 
2 ie se ew £8) (64 6) 6: a) Je: 


AANNNNN Semen 


ODNSIMO Daa 
Mmrncsceoonamare 
eoervreereeeeeee 
AD 3Z00~—e ee 


ON UOMG ~ RONG 
CCeNR arr oNns 
eoeeteeee et eee 


OK NNMMaASSt TS 


oscooceoccooo 
RNNANNNNAN 
oot eevee eee 


MAMAMAM AMA AM Moy 


MAMA MA AM amas 
et et tt tt ot 
eer ee ee eevee 


PAAPRAAARAAAN 


AOPDOM=DOMN 
DAP rsr mamas 
eererererereev ee 


Talololalelolololojo) 


DTFRANAP MAMAS 
OM INTFANODOS © 
eseeeereese 


OM ANNMTOO OF 


mee He 
wloleloelololelolelajo} 
eoeereeeveeee 


lelalolelslolslalolols) 


WADIA FPP 
SOAMe AA MGMmm 
e: 2 19 18 @ oe (@ 6 6.6 8 


ocoeoocoecoscoe 


NANA DPLPOOMAW 
ITPNOO-DODATO 
Se 2 e@ 8 8 © 66 @ 6 6 
zaroocnmocecc- 


OMOr Damme 
ONMNVOoOreomnanans 
oeeeerteetrtee 


COOKS =SNNMGmS 


HANA DOLAAE OMAN 
QOTOKHNAT DDH Sf 
eee ee eee ee om 
aq00s000000=C + 

7 

co 
OLFTUIMAMO em tt It 
OFC NmAsrOe Nore 


est eeeeee 2 CO 
COOCOOCOCCoO°O=$a=2 
$8 


SCAANINANFONOO YO 
NINNNNNN NNN 
eeeeeeeeeeve ~ 
lolslojalolojolainiclojal,. 
ou 
oo 
APN a AAMN eM & 
ITO OVOEEOTTOrO 


eee ee tee ee eDil 
lelelalelolololoiniol nae) 
an 
a- 


COBDWMO-EMANOA~ Trot 
OSANANMTNOADOE OC 
eoeoevreeereaeeee WY 
lol ole lolelololele Lalo) in 

On 
CONENYDMOOON e—« 
OnOsdTroanqgyonoe 
esos eeereee vet @ 
990 4 SN NOPE OD 

ww 

t-_ 
aNMNOO-FOTE DMO 
DAW OOP DF aeMOOO 


soeeereree eevee =] 
oCo0070700%-=——— 
a 
“” 


AMADOR WWNDAMPNOM 
WR ONDONSEOT om 


oeereereere ee Ov @ 
CRI NNMPPNOF 19 
awh 

“== 

eee e © & 8 oe Ome 
CaN OO™ORFOOQQ 
OU 


CA XCPW XCPT 


CN 


CLTD CODID 


CLIP CLUT CLB COL CNwP CNTP 


CLBwW 


PENS eR NUOW 
ANAK RENO IAN 
Oe) 8 67a © Jee) 16S 0 


PLTAANAMMOANN 


FROSILMOONT AA 
FOC noocrasquam 
eee eer eee eo 
MIANOD OER 


ONOM AE AOHNN 
CORMFaetrrony 
60) e@ 8 ‘e 2's e080 © @ 


OR ANN SIMMS PS 


ololalalalelolololoyo) 
NNNNGINNN NN 
OO 458) 16 6/19 10 BD) 6 6. 


AMAA AMM MAME 


ANTM AM AMAA AM 
et et eg et ot i tt 
eeoeeeevertree*e* 


FATAAPARAAAH 


NNAKAPE ANOS 
AUN a rr 
eoeeteeeeetee ® 


lelolelplalelelelolola) 


NONASHNO TO 
AMTNAONNOOE 
eeoeeereeweetee @ 


ORANG PNOOF 


et et tt tet 
lolelolelojlolelelajleloa) 
ereev eevee ee % @ 


on00e 71000900 


DWODNMNNMSPmS 
rere 
eevee rtee rave 


scooosocoeseo 


FAND 2H LFOUoN 
NOSE DORON 
eeoeveeevre 8 @ 
fsleleipialelal l.lal 


QMOromor4nyro 
OmoOoOrenamam 
eeeteoeevee eee 


COOK AANNMO SS 


SOSH OE OM MIA 
DDO Ra PHOy 


eee eee ee eo orm 
esa 300970700=0 
ee 
foto} 


OILAMAMMAMOME™ 1H 
DOC ANASAE Dele 


eestor eee & CM 
lolalololelololololol. [a> 4 
Vo 

a 


CH DOADOIOAH—N UO 
NNNNNN AMMA oD 
eee r ee eee ete 
lelelel=lolnial@lelelolol,:| 
ow 
oo 
CWA SFOONOLT ee 
WOE SAHROCaHNEO 
a ee) 8) 6) ase. « 2) 8Co 8 
lolelelelaloloh lalla: 40) 
an 
iol 
CONAOCEANCETOET 
ODOHeINIAMS FAN OS OC 
eeovoeveeveevee eevee WY 
B97090900909009090 in 


On 
ONEANDMOTA tae 
ONOCOLlTNONE NO 

eereeever eevee @ 

Ihe 
NIE FORLENTOQO 
WINWOOD DONFEUO 

sere everee eee WL 
O00090090-——--<4 

-_ 

aa) 
WWNO Dae ~eaHDPProm 
OMAP aE PTNUTMHO OM 
o*ere ee ees ae 8 @ 
ORS NMOMTANOF HO 

—t 

zt 


—“ 


Table (V-4). (Continued) 


81 


iP Ss OP Tage amma Pe (Gln Se ak Os (Rca aa anmcee ) " ae OAT a lane Sa Pa mae.) Coane ate he aaa 
i ' ' ' \ i 
i i ' i i i 
Seco | i i i ' i 
' ' ' ' ‘ ' ' 
q oxree LIE SOI pt bees (Soe as ae ae) (a A (eee Tassos mag (eae cae Rass r ff aa 
i i i \ i i i i 
t i { i ' i i i 
i i | i 
i \ i ! ! | i 
farina | (Rc r es ee a a Sal 
‘ i i i | ' 
‘ 1 f ' i \ | i 
i i 1 ' ' \ ' 1 
i i i \ i i i i 
(eee eal (ea aaa Cae ie ee haa it 
i i i i ' i i i 
i i i i i i i i i 
i i i i \ i i i i 
i i i i i i i i i 
Rate ss Sr = fae el Baer ans ce dle i Ome aes) Ate 7 See eee if 
\ i i i { i i 
i i i 
‘ ' ' ' ' ! ' 
i i i i i \ i 


eee ee: 
i 
1 
! 
i 
! 
4 
i 
i} 
' 
1 
I 


on eeteetcteriate bed 
' 

i 

! 

i 


! 
1 
i 
! 


t 
1 
1 
1 


SS OS 


a 


Output data plot for Example V-B 
82 


Figure (V-4). 
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Table (V¥V-5). Output variables of LAERO1 
Angle of attack 
Total coefficient of litt 
Total coefficient of drag 
Wing panel coefficient of lift 
Additional [ift on body due to wing 
Tati apane! coefkticient of Litt 
Additional Fitt on body due to tail 
Body alone lift coefficient 
Induced drag coefficient 
Wing panel normal force coefficient 
Tail panel normal force coefficient 
L2f£e cocttictent due to tail deflection 
Drag coefficient due to tail deflection 
Total normal force cosfficisnat 
Total axial force coefficient 
Wing center of pressure 
Tail center of pressure 
Total missile center of pressure 
Total pitching moment about C.G. 


Zero lift drag coefficiant of wing-body- 
tail Scomb iadt ron 


Miscellaneous zero lift drag coefficient 
zero PEt idceag “coefficient ‘of tail 

zero lift drag coefficient of wing 

Zero lift drag coefficient of body alone 
Drag coefficient of body protrusions 
Drag coefficient of engine inlet 

Drag coefficient of boattail region 
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Sider Ts 


VI. CONCLUSIONS AND RECOMMENDATIONS 


There are many topics which may be the subjects of fol- 
low on work contained within this thesis. Although the four 
programs have been installed on the IBM 370 computer system, 
these four alone do not fully satisfy the original goal of 
this work: Provide a computer supplement to the Tactical 
Missile Conceptual Design handbook. Numerous additional 
focal algorithms are utilized in the design handbook which 
deserve the attention of a programmer. Of immediate inter- 
est are the areas concerning radar or infrared guidance 
systems, baseline configuration modelling and weight distri- 
bution, and initial control and lifting surface design. 

Each of these topics can be programmed to provide missile 
design students interactive learning tools when coupled with 
the design handbook. 

The most urgent follow on work to this thesis is the 
restoration of the program LAERO1 to a reliable, useful pro- 
gram. The program was modified and set up on the IBM 370 
computer system during the period immediately following the 
System's installation at the Naval Postgraduate School. As 
could be expected, the computer suffered many and varied 
grcwing pains in its early life. As a result of this, or of 
the human manipulation expanding the capability of the pro- 
gram, the effectiveness of LAERO1 was substantially reduced. 

Work involving the other three programs would involve 
Simply expanding their capabilities. The trajectory models 
program, LPATH, presently considers only two guidance laws: 
line-of-sight guidance and proportional navigation guidance. 
Other guidance laws which can easily be included in the pro- 
gtam include pursuit, beam rider, and combinations of diffe- 
Tent laws. It might also prove useful to be able to 
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simulate the entire missile trajectory but still only output 
the terminal phase of the encounter. Another option would be 
tc provide the target with a controlled trajectory instead 
cf the constant acceleration condition now imposed. 

The propulsion sizing program, LPROP, should have the 
various nozzle options incorporated into the program so that 
it isn't necessary to manually juggle the program output. 
Other booster-sustainer grain configurations could be 
explored, such as the booster grain being cast within the 
ccre of the sustainer, or even a motor with only a single 
grain. Another suggestion for the convenience of the pro- 
gram users is to institute a shopping list of availlable 
prepellants and their characteristics into the progran. 

A similar list of available explosives and case materi- 
als could be put into the warhead sizing program, LBOMB. 
These shopping lists would provide ready access to reference 
information and, at the same time, decrease the number of 
data values to be manually input into the computer. Since 
the current program is limited to cylindrical warheads, an 
area of expansion would be the flexibility of warhead 
styles, such as curved, shaped charge, continuous rod, etc. 

At present, the programs are somewhat hindered by the 
mechanics involved in producing the printouts and the plots. 
Due to the results of tailoring a program to be interactive, 
often it is required to completely exit a program before 
Output can te received. Subsequent design iterations 
Teguire re-entry into the program, which produces a certain 
justifyable annoyance to the user. Additionally, the compu- 
ter center has instituted a new policy of cancelling any 
jobs with duplicate job names, which can be severely irri- 
tating and cumbersome to the persons running the plot rout- 
ines contained within this thesis. However, the computer 
center has developed procedures which have the potential to 
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alleviate both of these problems. According to Volume 13, 
Number 4 of the Computer Center Newsletter, CMS commands can 
now be invoked from within a FORTRAN program. The print and 
plot operations presently contained within executive rout- 
ines can now be placed directly within the source programs. 
This will remove all current restrictions placed om the nun- 
bers of printouts received per session and will label each 
plot with the user's job name and not the project's job 
namre. 
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APPENDIX A. TRAJECTORY MODELS PROGRAM FLOWCHART 


T, TSP, TAC, TAL, MSP,IR 


Input data: TITLE,DEL,NAV,LOS 
IR f, f 


Playback input data to terminal 
for veELincation 


TTYI=LOS+TAL 
MTH=LOS+MAL 


LOS=DETORA(LOS 
TAL=DETORA TAT 
TTH=DETORA(TTH) 
MAL=DETORA (MAL) 
MT H=DETORA (MTH) 
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A2=arsin (sin (A1) *TSP*IRM/(4SP*IRT) ) 
MAL=A2 
MTH=MAL+LOS 
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(D5+D7+D8) /D9 
T1*+DEL*T3 
P1=P1+D EL*P3 
R2=R2+DEL*R4 
A2=A2+DEL*A4 


F2 


A1+#DEL*A3 
R1=R1+D EL*R3 
T2+DEL*T4 
A(I) =MSP*T4 


T1= 
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A4= 
Al 
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+DEL*MSP*cos 
+DEL*MSP¥sin 
+DEL*TSP*cos 
YT +DEL*TSP*sin 
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A2=arsin (sin (A1) *7 
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A2 
MTI=MAL+LOS 

=0 

=0 


ERT<-ERE 


MAL 
R (I) = 


arsin (sin (A1) *TSP/msp) 
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*Sin 
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=XM (I-1) #*MSP*cos (T2) *DEL 
YM{I) =¥M (1-1) #+MSP¥*sin(T2) *DEL 
XT (I) =XT (1-1) +TSP*cos(T1) *DEL 
YT (I) =Y¥T (I-1) #TSP*sin(T1) *DEL 


[ _R2=7SP*cos (A1) -MSP*sin(A2)_ 


noe a Cesta VAESBATAC# COs (41) 78 (T-1) eae 
D9=A (I-1) ® (2*R2+NAV*MSP¥coS (AZ) ) /R (I-1) 


[ a3=8-p9 | 


| D2=(TSP*sin (Al) -4SP*sin (A2)) /R(I-1) | 
D6=TAC/TSP 


D7=NAV*D2 
P2=D6-D2 
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R(1L) =R (I-1) +D EL*R2 
A({I) =A (I-1) +DEL*A3 
P1=P1+DEL*D2 
T1=f1+DEL*D6 
T2=[2+DEL*D7 
A1=A1+tDEL*P2 
A2=A2+DEL*P3 


Call SCREEN 
Do you want to run the 
same problem again? 
DE (LOS) 
DE (EAL 


Do you want -o run \ = 
a new problen? 
[ End | 
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APPENDIX B. TRAJECTORY MODELS PROGRAM LISTING 


Following the next page is the program listing for the 
Trajectory Models program. It consists of three main divi- 
sions; the executiv2 routines, the FORTRAN IV computational 
program, and the FORTRAN IV plotting program. The 2xecutive 
routines are used to establish required file definitions, 
initiate operation of the computational program, supervise 
the transfer of data to the plotting program, and provide 
operational information to the program user at appropriate 
times. 

The computational program, LPATH FORTRAN, consists of 
four subprogram divisions. The MAIN program accepts the 
input data and performs th2 calculations for the line of 
sight and the proportional navigation problems. Subroutine 
SCREEN is used to prompt the user to clear the terminal 
screen for proper positioning of the output. Subroutine 
SORT determines the largest missile acceleration value and 
the value ranges of the X and Y position coordinates for 
plotting reference. Finally, subroutine OUTPUT formats the 
Calculated data and provides it to the user, the printer 
file, and the plot program data file. 

The plot progran, PATHPLOI FORTRAN, is structured for 
direct submission to the MVS batch reader from the 3278 ter- 
Minal. No cards need to be punched or read. Attention must 
be given to the correct JCL accounting data in the first 
line of the program; thos2 shown in the listing are for 
illustrative purposes only. This program receives the data 
from the computational program and produces a singl2 geo- 
graphic chart of the encounter in the encounter plane. Mul- 
tiple problems, up to ning will overlay on the single 
chart. 
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APPENDIX C. WARHEAD DESIGN PROGRAM FLOWCHART 
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BU= (2*CTHS (CDEN/XDEN +1) *FRM 
* (L WH+ 2*RKILL*tan (B 
BL=5*CD EN* PKA¥ WT 
RW= ((-BU/BL) *alog(1-PD (I 


XW= (PI* (RW-CTH) **2) *LWH*X DEN 
CW=PI*(RW* ir i Cd 
HW=KW+CW 
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RMLB({1 Z§ one CODEN) iti 
HITS=FRAGS*A 
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APPENDIX D. WARHEAD DESIGN PROGRAM LISTING 


This program has two major sections; the executive rout- 
ine and the FORTRAN IV computational program. The executive 
routine establishes the required file definitions and initi- 
ates operation of the computational progran. 

The computational program, LBOMB FORTRAN, consists of 
four subprogram divisions. The MAIN program accepts the 
input data, calculates the atmospheric characteristics, and 
formats and displays the output to the user and sends it to 
the printer file. Subroutine VEL calculates the initial 
velocity required to propel a given mass a specified dis- 
tance through the atmosphere with a particular residual vel- 
ocity remaining. It also jietermines the ballistic limit 
velocities for the situation. Subroutine SIZER siz2s the 
warhead for a given Pd value. It also produces the charge- 
to-mass ratio, the number of fragments, the fragment size, 
and the average number of hits received by the targst. 
Subroutine SCREEN prompts the user to clear the terminal 
screen for proper positioning of the displayed data. 
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APPENDIX F. PROPULSION SIZINS PROGRAM LISTING 


Following this page is the program listing for the Pro- 
pulsion Sizing Program. It has two segments; the 2xecutive 
routine and the FORTRAN IV computational program. The exe- 
cutive routine establishes the required file definitions and 
initiates operation of the computational program. 

The computational program, LPROP FORTRAN, consists of 
six subprogram divisions. The MAIN program accepts the 
input data from the user and performs the guiding calcula- 
tions for the booster and sustainer motors. Subroutine 
PRATIO determines the ambi2nt pressures at the design alti- 
tudes. Subroutine PCRAT defines the optimum chamber pres- 
sure to ambient pressure ratio with respect to the case 
Material properties. Subroutine ARAT solves for the area 
ratio of the nozzle and tries to size the nozzle to the mis- 
Sile diameter by varying the chamber pressure, characteris- 
tic velocity, and thrust coefficient. Subroutine SCREEN is 
used to prompt the user to clear the terminal screen for 
proper positioning of the sutput. And subroutine OUTPUT 
formats the calculated solutions and provides them to the 


user and to the printer file, if so directed by the user. 
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APPENDIX G. AERODYNAMIC COEFFICIENTS PROGRAM LISTING 


This program is divided into three major subdivisions; 
the executive routines, the FORTRAN [I/IV computational pro- 
gram, and the FORTRAN IV plotting program. The executive 
routines are used to establish the required file defini- 
tions, initiate operation 9f the computational program, 
supervise the transfer of iata to the plotting program, and 
provide operational information to the user at appropriate 
times. 

The computational program, LAERO1 FORTRAN, consists of 
five subprogram divisions. The MAIN program accepts the 
input data, conducts calc lations not done by other subrout- 
ines, formats the calculated data and provides the sutput 
data to the user, the printer file, and the plot program 
data file. Subroutine GEOSUB calculates the missile wetted 
area and the Reynolds number per foot based on the flight 
altitude. Subroutine CLASJB calculates the aerodynamic sur- 
face lift-curve slopes. Subroutine CATSUB calculates center 
of pressure locations, cross-flow drag coefficients, and 
interference factors. Subroutine SCREEN prompts the user to 
Clear the terminal screen for proper positioning of the out- 
put. 

The plot program, AEROPLOT FORTRAN, is structured for 
direct submission to the MVS batch reader from the 3278 ter- 
Minal. No cards need to ba panched or read. However, 
attention must be given to the correct JCL accounting data 
in the first line of the program; those shown in the list- 
ing are for illustrative purposes only. This program 
receives the data from the computational program and pro- 
duces a group of six charts for each Mach number entered. 
The plots represent the relationships o£ Cl to alpha, Cn, 
Cn, Cd, Ca, and Cdi. The program can produce up t> 24 sets 
ez plots for a single run. 
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FILE: LAERO1 EXEC A NAVAL POSTGRADUATE SCHOOL 


FILEDEF 08 DISK LAERO1 OUPPUL AD (RECFM VA BLOCK 133 PERM 
Beer OF DISK LAERO1 PLOT AO (RECFM VA BLOCK 80 PERM 
YOU WILL HAVE THE OPPORTONITY TO OBTAIN BOTH A HARDCOPY 
PRINTOUT AND A SET OF PLOTS FOR ONE SET OF INPUT PARAMETERS 
EACH TIME YOU ENTER THIS PROGRAM. THE PROGRAM BE 
RE-RUN CONTINUOUSLY AND YOO WILL HAVE THE OPTION [90 CHANGE 
INPUT PARAMETERS FOR EACH SUCCESSIVE RUN, BUT YOU CAN 
OBTAIN THE PRINTOUT AND PLOTS PERTAINING TO THE LAST RON 
ONLY. IF ADDITIONAL OUTPIT IS REQUIRED, RE-ENTER THE 
PROGRAM. 
SEND 
LOAD LAERO1 
START 
& BEGTYPE 
TO OBTAIN THE HARDCOPY PRINTOUT OF THE DATA TABLES, TYPE 
AND ENTER: 
LAERO1PR 
Eateee THE VERSATEC PLOT OF THE TABULAR DATA, IYPE AND 
LAERO1PL 
& END 
PIGe:« LABROIUPR» .EXEC A NAVAL POSTGRADUATE SCHOOL 
PRINT LAERO1 OUTPUT (LINEZSUN 70 
&BEGTYPE 
THE OUTPUT WILL BE PRINTED ON THE VM PRINTER IN ROOM 140 
AND WILL BE IDENTIFIED WITH YOUR USER NUMBER AND LAST NAME. 
so a eae REQUIRES “SEVERAL SINUTES TO OBTAIN THE PRINTOUT. 
FILE: LAEROIPL EXEC A NAVAL POSTGRADUATE SCHIOL 
COPY LAERO1 PLOT A LAERO PLOTDATA A PLOT LAERO1 A 
EXEC SUBMIT PLOT LAERO1 A 
ERASE PLOT LAERO1 A 
SBEGTYPE 
THE PLOT WILL BE DRAWN IN THE COMPUTER ROOM AND PLACED OVER 
fae ALPHABETIZED OUTPUT FILE IN ROOM 140. IT WILL. BE 
IDENTIFIED BY "AEROFILE" AND USUALLY REQUIRES MANY MINUTES 
pomp DAYS) TO OBTAIN. 
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